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ABSTRACT 


-  ; 

Adaptive  antennas  provide  an  important  means  of  enhancing  signal-to-noise 
ratio  in  the  adverse  electromagnetic  environments  that  sometimes  arise  due  to 
jamming  or  interference.  For  that  reason,  adaptive  antennas  are  increasingly 
finding  application  in  high  performance  radar  and  communications  systems.  In 
many  of  these  systems  the  adaptive  subsystem  must  exhibit  robust  performance 
in  the  face  of  multipath  or  rapidly  changing  interference. 

Unfortunately,  the  adaptive  beamformers  now  in  use  do  not  perform  well  in 
certain  environments.  It  has  been  known  for  some  time  that  correlated-signal 
conditions  (e.g.,  due  to  multipath)  can  lead  to  partial  or  total  cancellation  of  the 

i 

desired  signal  within  *:n  adaptive  beamformer.  More  recently  it  has  been 
recognized  that  signal  cancellation  can  arise- during  high-speed  adaptation  even 

though  1)  the  desired  signal  and  interfering  ^signals  are  uncorrelated  and  2)  the 

/  • 

look-direction  response  is  rigidly  controlled.  The  initial  sections  of  this  report 
examine  the  signal- cancellation  mechanism,  with  emphasis  on  the  more  difficult 
uncorrelated-signal  case.  Simulation  results  are  presented  that  illustrate  the 

i  . 

cancellation  effect,  and  an  analysis  is  given  for  a  simple  environment  consisting  of 
the  desired  signal  and  one  jammer.  '  " 

The  remaining  sections  of  the  report  describe  a  technique  for  avoiding  signal 
cancellation.  The  adaptation  problem  is  reformulated  to  permit  jammer  nulling 
and  signal  recovery  under  conditions  that  ordinarily  result  in  signal  cancellation. 
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The  reformulated  problem  suggests  a  new  adaptive  structure  comprised  of  a 
master  beamformer  in  which  adaptation  is  conducted  in  a  synthetic,  signal-free 
environment  and  a  slaved  beamformer  in  which  the  actual  array-element  signals 
are  processed.  The  new  structure,  which  is  termed  a  composite  beamformer,  is 
compared  to  Frost’s  hard-constrained  adaptive  beamformer.  Particular  attention 
is  given  to  relative  performance  in  cancellation  environments,  to  conditions  under 
which  optimal  behavior  is  approached,  to  convergence  speed,  and  to  weight 
behavior  near  convergence. 
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L  INTRODUCTION 


The  past  two  decades  have  witnessed  the  emergence  of  the  adaptive  antenna 
as  an  important  element  in  countermeasures-resistant  radar  and  communications 
systems.  The  application  that  has  inspired  most  of  the  development  effort  is 
jammer  nulling;  adaptive  antennas  can  eliminate  a  substantial  fraction  of  the 
incident  jammer  power  while  admitting  some  desired  signal  (or  signals)  with 
useful  gain.  A  key  aspect  of  this  selective  treatment  of  jammers  and  desired 
signals  is  that  it  is  accomplished  without  extensive  a  priori  knowledge  of  the 
jammer  environment.  Instead,  the  locations  of  jammers  are  “learned”  by  some 
adaptive  algorithm,  and  nulls  are  automatically  steered  to  and  maintained  on  the 
jammers. 

One  of  the  fundamental  issues  in  adaptive  antennas  has  been  signal 
preservation.  It  has  been  clear  from  the  beginning  that  the  nulling  behavior 
which  is  of  such  great  value  in  eliminating  jammers  is  a  two-edged  sword  that 
also  threatens  the  desired  signal.  Proper  control  is  absolutely  essential  to  avoid 
signal  loss  in  the  adaptive  beamformer  (ABF).  The  original  sidelobe-canceller 
scheme  of  Howells  and  Applebaum  (1,2]  exploited  the  differing  signal-to-jammer 
ratios  in  a  directive  primary  antenna  and  an  omnidirectional  auxiliary  antenna  to 
avoid  seriously  attenuating  desired  radar  signals.  Widrow  [3],  operating  with 
somewhat  different  signal  and  antenna  assumptions,  took  the  tack  of  introducing 
pilot  signals  to  control  beamformer  response  in  specified  look  directions.  Griffiths 
[4]  devised  a  different  soft-constraint  technique  that  involved  statistical 


characterization  of  the  desired  signal'-  rather  than  the  actual  surrogate  signals 
required  by  the  pilot-signal  scheme.  Frost  [5,6]  developed  a  constrained  least- 
mean-square  (LMS)  algorithm  that  assured  exact  conformance  with  some 
prespecified  look-direction  response.  More  recently,  Griffiths  and  Jim  [7,8] 
contributed  a  structure  called  the  “generalized  sidelobe  canceller,”  which 
provided  an  alternative  method  of  realizing  hard  constraints.  Some  further 
development  and  generalization  of  soft-constraint  methods  has  also  taken  place 
within  the  past  few  years;  Chestek  [0]  brought  together  much  of  the  earlier  work 
on  soft-constrained  methods  by  combining  soft  linear  constraints  with  a  mean- 
square-error  criterion  in  his  soft-constrained  LMS  algorithm. 

A  unifying  factor  in  the  work  just  outlined  is  the  focus  on  the  response  of 
the  antenna  to  desired  signals.  This  approach  involves  a  tacit  assumption  that 
the  problem  decomposes  in  a  tractable  way,  i.e.,  that  desired  signals  and  jammers 
can  be  treated  independently.  The  assumption  can  be  justified  in  many  cases, 
but,  as  will  be  demonstrated,  there  are  some  rather  simple  signal/jammer 
scenarios  where  difficulty  arises.  When  a  failure  occurs,  the  usual  constraint 
methods  do  not  adequately  protect  the  desired  signal.  Instead,  the  signal  is 
partially  or  totally  destroyed  in  the  adaptive  beamformer  by  residual  jammer 
signals.  This  phenomenon  is  termed  signal  cancellation. 


The  remainder  of  this  report  is  devoted  to  a  further  discussion  of  signal 
cancellation  and  to  the  description  of  a  method  for  avoiding  it.  The  next  chapter 


provides  insight  into  the  signal  cancellation  phenomenon  by  explaining  how  it 
was  recognized  and  by  showing  the  effect  in  two  simple  cases.  Chapter  m 
introduces  a  composite  beamformer  that  defeats  signal  cancellation  and 
demonstrates  the  performance  improvement  that  it  provides.  Chapter  IV 
elaborates  upon  the  characteristics  of  the  composite  beamformer  by  describing 
convergence  behavior,  illustrating  performance  under  additional  signal  conditions, 
and  examining  the  issues  of  convergence  rate  and  weight  behavior  near 
convergence.  Finally,  Chapter  V  summarizes  the  status  of  the  work  on  signal 
cancellation  and  lists  some  topics  that  deserve  further  research. 


D.  SIGNAL  CANCELLATION  IN  HARD-CONSTRAINED 


BEAMFORMERS 

Jamming  or  interference  generally  involves  the  addition  of  spurious 
components  to  the  desired  or  target  signal.  These  spurious  components  can 
decrease  intelligibility  or  cause  decoding  errors,  sometimes  to  the  extent  of 
denying  communication  or  surveillance  altogether.  Adaptive  antennas  are  useful 
in  suppressing  the  “additive”  jamming  effect,  but  can  be  responsible  for 
introducing  a  “subtractive”  jamming  effect.  Subtractive  jamming  amounts  to  the 
cancellation  of  desired  signal  components  and,  like  ordinary  jamming,  can  disrupt 
communications.  The  purpose  of  this  chapter  is  to  provide  insight  into  the 
nature  of  this  signal  cancellation  phenomenon. 

The  cancellation  phenomenon  will  be  examined  in  the  context  of  adaptive 
beamformers  that  use  versions  of  the  LMS  algorithm.  It  will  be  seen,  however, 
that  cancellation  arises  not  from  any  special  properties  of  the  adaptation  method, 
but  due  to  nonzero  correlation  between  the  signal  and  jamming  waveforms 
during  the  data  window  being  considered  by  the  adaptive  algorithm.  Similar 
cancellation  effects  will  take  place  with  any  least-squares  algorithm  that  might  be 
used  to  adjust  the  beamformer  weights  at  a  comparable  rate. 

The  chapter  has  been  divided  into  several  sections.  The  first  section 
furnishes  a  bref  historical  perspective  on  signal  cancellation.  The  second  section 
pro”' J  »  a  /eview  of  the  properties  of  Frost’s  hard-constrained  algorithm  and 
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establishes  notation  that  will  be  useful  throughout  the  report.  The  third  section 
discusses  the  effects  of  signal  cancellation  on  a  wideband  signal  and  illustrates  the 
destruction  of  signal  components  that  can  occur.  The  final  two  sections  consider 
cancellation  in  the  narrowband  case.  Appendix  A  contains  a  supporting  analysis 
for  the  final  section. 

A.  Background 

It  has  been  realized  for  some  time  that  signal  cancellation  can  occur  even  in 
hard-constrained  beamformers  when  correlation  exists  between  the  desired  signal 
and  some  other  signal  impinging  upon  the  array.  (See,  for  example,  [6]  by  Frost.) 
This  was  a  disquieting  fact,  because  this  is  precisely  the  condition  that  can  arise 
naturally  as  a  result  of  multipath  or  can  be  induced  artificially  through  repeater 
jamming.  It  appears,  however,  that  this  troublesome  case  was  more  or  less  set 
aside*  while  work  continued  on  other  aspects  of  adaptive  antennas. 

As  adaptive-antenna  work  went  forward  in  other  areas  during  the  mid¬ 
seventies,  developments  relevant  to  signal  cancellation  were  taking  place  in 
another  branch  of  adaptive  systems.  Experimental  work  in  adaptive  noise 
cancellation  had  disclosed  filtering  phenomena  that  fell  outside  the  purview  of 
Wiener  filter  theory.  Glover  [10,11,12]  analyzed  this  non-Wiener  behavior  and 

*Some  work  directed  specifically  towud  adaptive  elimination  of  multipath  interference  at  HF  waa  described  in 
November  1981  by  Hansen  and  Lonfhlin  |19].  Their  method  makes  use  of  a  modulated  pilot  signal  that  is  added  to  the 
communication  signal  at  the  transmitter.  Pilot  bandwidth  mutt  be  adequate  to  provide  discrimination  between  the  vari¬ 
ous  modes  arriving  at  the  receiver.  This  approach  is  limited  to  those  sitnations  where  the  transmitter  it  accessible  and 
bandwidth  allocations  are  not  restrictive. 
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showed  that,  for  the  case  involving  sinusoidal  reference  inputs,  the  adaptive  noise 
canceller  can  be  regarded  as  a  stable  notch  filter.  Shensa  [13,14]  was  later 
successful  in  treating  a  somewhat  more  general  class  of  input  signals. 

Widrow  was  able  to  extrapolate  from  the  work  being  done  on  adaptive  noise 
cancellers  and  draw  some  disturbing  new  conclusions  concerning  signal 
cancellation  in  adaptive  beamformers.  Widrow  realized  that  correlation  (in  the 
usual  long-term  sense)  was  a  sufficient,  but  not  necessary,  condition  for  the 
destructive  effects  of  signal  cancellation.  If  high  adaptation  speeds  are  employed, 
signal  cancellation  can  be  induced  by  a  broad  class  of  jammer  signals,  and  close 
replication  of  the  desired  signal  (e.g.,  through  repeater  jamming)  is  not 
necessarily  required.  Since  high  adaptation  rates  are  being  sought  so  that  good 
jammer  nulling  performance  can  be  obtained  in  dynamic  countermeasures 
environments,  the  conditions  that  support  signal  cancellation  may  be  routinely 
present.  Recognition  of  this  state  of  affairs  added  considerable  incentive  to  study 
signal  cancellation.  Simulation  work  was  begun  at  Stanford  in  1080  to  confirm 
the  effects  predicted  by  Widrow  and  to  support  the  search  for  methods  of 
avoiding  signal  cancellation.  Portions  of  this  work  have  been  reported  in 
References  15  and  16.  Other  workers  [17,18]  have  observed  signal  cancellation  in 
actual  adaptive-antenna  systems. 


B.  Frost's  Constrained  LMS  Algorithm 

The  application  of  his  constrained  LMS  algorithm  to  adaptive  arrays  has 
been  described  in  some  detail  by  Frost  [5,6].  Some  of  the  key  results  will  be 
restated  in  this  section  to  provide  the  background  and  notation  for  discussions 
that  follow. 

The  adaptive  array  problem  to  be  considered  is  pictured  in  Figure  1.  A 
uniform  linear  array  of  K  elements  is  connected  to  an  adaptive  beamformer 
that  consists  of  K  tapped-delay-line  filters,  each  with  J  taps  and  J 
adjustable  weights.  The  sampled  tap  voltages  and  the  weights  are  indexed  in  the 
columnar  scheme  shown  in  the  figure  and  are  identified  in  vector  notation  as 

xTU)  4  !*,</),  *s0), .  **,(/)]  (2-1) 

WTU)  4  |w,(»,  »j|j) .  wH|j)|  ,  (2-2) 

where  j  is  the  sample  number.  The  filter  outputs  are  summed  to  form  the 
beamformer  output: 

vU)  =  WTU)X(i)  =  XTUWU )  .  (2-3) 

In  addition  to  being  the  useful  output  signal  from  the  system,  the 
beamformer  output  serves  as  the  error  signal  that  is  fed  back  for  use  in  the 
weight  adaptation  process.  The  goal  of  adaptation  is  to  minimize  the  output 
contribution  of  noise  sources  (such  as  the  jammer  indicated  in  the  figure),  subject 


Equivalent  Look-direction  Filter 


Diagram  of  the  Frost  Adaptive  Beamfo rmer 


to  a  set  of  constraints  that  provide  rigid  control  over  the  array  response  in  the 


direction  from  which  desired  signals  are  expected. 

It  is  assumed  that  the  desired  signal  and  the  interference  are  uncorrelated. 
The  incident  signal  vector  may  be  decomposed  according  to 


xu)  =  su)  +  mi) 


(2-4) 


by  letting  S(j)  represent  the  vector  of  desired  signals  at  the  beamformer  taps 
and  N(j)  represent  the  aggregate  contribution  of  the  various  noise  sources.  If 
the  autocorrelations  of  the  quantities  of  Eq.  2-4  are  defined  by 


Rxx  4  EWi)  XTU)] 

Kss  4  £|S(»  Sr<»] 

Rm  4  *W)  NTU$  , 


(2-5) 

(2-6) 

(2-7) 


then  the  assumed  lack  of  correlation  between  the  desired  signal  and  the  noise  is 
expressed  by 

Rxx  =  Rss  +  Rnn  (2*8) 

This  correlation  condition  must  hold  if  the  look-direction  constraints  are  to  afford 
protection  for  the  desired  signal. 

The  constraint  method  depends  upon  knowing  the  arrival  direction  of  the 
desired  signal.  The  array  is  assumed  to  be  steered  (either  mechanically  or 
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electrically)  toward  the  known  look  direction.  Under  this  assumption,  the  desired 
signal  components  are  identical  at  the  inputs  to  the  beamformer  filters  and  along 
any  column  of  taps  in  the  beamformer.  An  equivalent  look-direction  processor 
can  therefore  be  formed  by  summing  the  weights  in  each  beamformer  column  to 
arrive  at  the  J-tap  filter  shown  at  the  top  of  Figure  1.  The  equivalent  look- 
direction  filter  provides  the  key  to  an  algorithm  that  minimizes  noise 
contributions  at  the  output  without  influencing  the  response  to  the  desired  signal. 
Specifically,  the  weights  may  be  freely  adjusted  to  minimize  the  noise  output 
provided  the  column  sums  in  the  beamformer  remain  equal  to  the  preselected 
weights  in  the  equivalent  look-direction  processor.  In  order  to  state  the 
minimization  problem  succinctly,  it  is  helpful  to  express  the  look-direction 
impulse  response  in  vector  form  and  to  define  a  constraint  matrix.  The  look- 
direction-response  vector  is  formed  from  the  weights  in  the  equivalent  look- 
direction  filter: 

l  4  l/i,  ft .  //]  .  (2-9) 


rggm'JCTE  «**■!■  >  J-"  rr  *■.■  ■r.r.".*.  tj  .■•  v-  .■•  r--*  v  •  •  \-  v  ■ 
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The  constraint  matrix  C  is  composed  of  columns  of  the  form 


c,  = 


|  1**  group  of  K  elements 

|  (t— 1)**  group  of  K  elements 
|  ingroup  of  K  elements 


(»  +  1)**  group  of  K  elements 


) 


Jth  group  of  K  elements 


That  is, 


C  ^  [cj,  Cj,  ...»  c»>  •••»  C/1 


(2-10) 


(2-11) 


The  matrix  C  is  a  KJ  X  J  matrix  that  conveys  the  beamformer  structural 
information  needed  to  implement  the  constraints. 

With  the  notation  developed  to  this  point,  the  constrained-LMS  problem 


may  be  stated: 
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Minimize  E[y2(j))  =  E[WTX[j)XT(j)\V]  =  WTRxxW  (2-12) 
subject  to  CTW  =  JT 

Frost  developed  both  closed-form  and  adaptive  solutions  for  the 
constrained-LMS  problem.  The  method  of  Lagrange  multipliers  was  used  to 
arrive  at  the  optimal  solution  in  terms  of  the  signal  statistics  and  the  constraints: 

w  -  R&C\CtR&C\-'  l  .  (2-13) 

This  solution  assumes  that  the  signal  statistics  are  known  (or  have  been 
estimated)  and  imposes  a  heavy  computational  burden  if  the  matrix  inversions 
must  be  repeatedly  performed  to  track  a  dynamic  environment. 

The  adaptive  algorithm  generates  an  asymptotically  optimal  solution  to  the 
constrained-LMS  problem  and  requires  much  less  computation  per  weight-update 
cycle  than  the  closed-form  solution.  The  update  scheme  is  most  easily 
understood  from  the  individual  weight  recursions: 

w  i(/+ 1)  =  w,(;)  -  f*y{j)*\U)  -  k(/)  -  wtiMJ)}  +  ^ 

w*(;+  1)  =  wj r(j)  -  vy(j)xK(j)  -  £  (u>,(;)  -  /iy(;>,(»j )  +  ^ 

w#f+i(;+i)  =  ™ K+iti)  -  t*y(i)*K+iU)  -  E  k(»  -  #**(»*(/)!  + 


w2j <U+ 1)  =  w2A-(»  -  tty{j)x2K{j)  -  s  K(»  -  »y(i)xiU)\  +  ^ 

*kjU+  i)  =  ^kjU)  -mU)*kjU)  £  KO)  -  vy(jUiU)\  +  -77  ■ 

A  (/-1)£+  1  A 

(2-14) 

The  right-hand  side  of  each  recursion  involves  the  current  weight  value  plus 
several  update  terms.  The  first  update  term  depends  upon  the  current  error  and 
the  signal  at  the  weight  in  question.  The  adaptation  constant  ft  is  used  to  scale 
this  term;  ft  is  set*  to  yield  stable  system  operation  and  to  provide  an 
appropriate  adaptation  rate.  The  second  update  term  is  the  average  correction 
over  a  column  of  weights,  and  the  third  update  term  is  a  fractional  allocation  of 
the  appropriate  weight  in  the  equivalent  look-directional  filter.  Taken  together, 
the  three  update  terms  drive  the  weights  toward  the  values  that  minimize  output 
power  subject  to  the  specified  linear  constraints.  It  is  easy  to  verify  the  efficacy 
of  the  constraint  method  by  summing  the  updated  weights  over  any  column  of 
the  beamformer  and  noting  that  the  constraint  is  always  satisfied. 

The  constrained-LMS  algorithm  may  be  written  in  matrix  form  to  provide  a 
more  compact  description.  It  is  helpful  to  introduce  a  /^./-dimensional  initial 
weight  vector  F,  where 

*  Factor*  foeerai»(  tht  choice  of  ft  are  ditciwed  is  Reference*  8,  8,  and  10. 
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F  4  C(CTC)~xl  ,  (2-15) 

and  a  KJ  X  KJ  projection  operator  P,  where 

P  4  /  -  C(CTC)~lCT  ,  (2-16) 

With  this  added  notation,  the  constrained-LMS  algorithm  may  be  stated  as: 

W(0)  =  F 

W(j+ 1)  =  P{W(j)  -  M(j)X(j))  +  F  .  (2-17) 


It  will  also  be  useful  to  note  the  time  constants  exhibited  by  the  hard- 
constrained  algorithm.  Frost  showed  that  the  matrix  P  Rxx?  determines  the 
rate  of  convergence  of  the  mean  weight  vector  toward  the  optimal  weight  vector 
W*.  Specifically,  convergence  of  the  mean  weight  vector  along  the  ia 
eigenvector  of  P  PjcyP  occurs  with  a  time  constant  given  by 


T;  — 


ln{l-ii<Ti) 


(2-18) 


where  <r,-  is  the  eigenvalue  associated  with  the  itk  eigenvector  and  where 
updating  at  the  sample  rate  (as  opposed  to,  say,  every  tenth  sample)  is  assumed. 


C.  Signal  Cancellation  in  the  Wideband  Caae 

The  effects  of  signal  cancellation  in  a  hard-constrained  ABF  can  be  readily 
demonstrated  by  simulating*  the  simple  adaptive- antenna  problem  diagrammed 
in  Figure  2.  The  selected  adaptive  antenna  was  a  two-element  Frost  ABF  with 

*  Ai  overview  of  the  (imnlitioo  methods  used  is  this  reeetreh  ie  lives  is  Appendix  E. 


two  weights  per  element.  A  wideband  desired  signal  S  with  power  Ps  equal  to 
one  was  incident  from  broadside,  while  a  narrow-band  jammer  J  with  power  Pj 
equal  to  10  was  incident  from  an  angle  of  30®.  The  ABF  was  constrained  to 
provide  unity  gain  and  a  flat  frequency  response  in  the  look  direction  by  setting 
the  constraint-vector  elements  /  j  and  /2  to  1  and  0,  respectively.  ABF  weights 
were  updated  using  the  constrained  LMS  algorithm,  which  in  this  particular  case 
may  be  written  from  (2-14)  as 

wi(i+ !)  =  j  lw i(>)  - /*sr(i)*i(i)]  -  ^  [w2(»  -  py(;>2(;)l  +  -^L 

w2  (;+  1)  =  “•  [W2  (;)  -  py(j)x£j)\  -^[w,  (;)  -  py{j)xx{j) ]  +  -j- 

w3(y+ !)  =  (w3(;)  -  /*y(;>3(;)]  -  ~  [w4(;)  -  nyU)**U)\  +  -y 

w4(>+ !)  =  \  [w4(;)  -  1  -  “  [w3(;)  -  i*y{j)xz(j)]  +  -j-  . 

(2-19) 

In  these  equations  p  is  the  constant  that  controls  the  rate  of  adaptation;  p  was 
set  to  1  X  10-2  in  this  simulation. 

Figure  3  shows  some  measurements  made  after  a  sufficient  number  of 
adaptations  to  allow  transient  phenomena  to  die  away.  Part  a)  of  the  figure 
shows  the  antenna  pattern  at  the  jammer  frequency  of  0.25  fttmp  and  confirms 
that  a  null  has  been  placed  on  the  jammer.  Part  b)  is  a  plot  of  the  frequency 
response  in  the  jammer  direction;  this  plot  also  shows  a  deep  null  at  the  jammer 
frequency.  The  third  plot  simply  confirms  that  the  specified  response  in  the  look 
direction  has  been  attained.  Taken  together,  the  three  plots  seem  to  indicate 
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Postconvergence  Plots  for  the  Four-element  Frost  ABF 
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b)  Jammer-direction  Frequency  Response 

c)  Signal -direction  Frequency  Response 


that  all  is  well  and  that  the  intended  objectives  of  signal  recovery  and  jammer 
rejection  have  been  attained. 

The  favorable  picture  provided  by  the  antenna  plots  of  Figure  3  darkens 
considerably  when  the  desired  signal  and  the  ABF  output  are  compared.  Figure 
4  is  a  plot  of  ensemble- averaged  spectra  of  the  desired  signal,  the  jammer  signal, 
and  the  ABF  output.  It  is  clear  that,  despite  very  rigid  control  of  the  look- 
direction  response  and  good  jammer  rejection,  the  desired  signal  has  not  been 
recovered  intact.  Instead  the  signal  has  experienced  frequency-dependent 
distortion.  Signal  components  on  the  skirts  of  the  signal  remain  unaffected,  but 
components  in  the  vicinity  of  the  jammer  frequency  have  been  totally  or  partially 
destroyed  in  the  ABF. 

The  extent  of  signal  cancellation  for  a  given  antenna  structure  and  signal 
environment  is  governed  by  the  rate  of  adaptation,  i.e.,  by  the  adaptation 
parameter  /i.  The  simulation  just  described  was  repeated  for  several  different 
values  of  /j;  output  spectra  for  the  series  of  experiments  are  assembled  in  Figure 
5.  Parts  a  and  b  of  the  figure  involve  lower  adaptation  rates  than  the  rate  used 
to  produce  the  results  in  Figure  4  and  reveal  lesser  amounts  of  signal  destruction. 
In  fact,  cancellation  is  almost  undetectable  in  Part  a.  Part  c  repeats  the  results 
shown  in  Figure  4.  Part  d  involves  a  large  n  (a  value  near  the  instability  point) 
and  shows  the  near-total  signal  destruction  that  can  occur  in  an  extreme  case.  It 
is  apparent  from  Figure  5  that  cancellation  occurs  over  a  considerable  range  of 
adaptation  parameter  selections  and  that  the  destructive  effects  can  be  quite 
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Figure  5.  Output  Power  Spectra  of  the  Four-weight  Frost  ABF 

for  Various  Adaptation  Constants 


important  when  high  adaptation  rates  are  sought.  It  is  also  clear  that,  although 
cancellation  effects  may  be  made  negligible  by  choosing  a  sufficiently  small  /*, 
orders  of  magnitude  in  convergence  rate  must  be  sacrificed  to  maintain  signal 
quality.  At  these  low  adaptation  rates  the  ability  to  track  a  dynamic 
environment  will  be  greatly  reduced. 

D.  Signal  Cancellation  in  the  Narrowband  Case 

A  wideband  signal  was  useful  in  the  examples  of  the  previous  section 
because  the  broad  spectrum  clearly  showed  the  extent  of  cancellation  effects. 
This  section  will  concentrate  on  developing  a  more  detailed  picture  c*  the 
cancellation  phenomenon  by  considering  the  internal  operation  of  a  beamformer 
while  cancellation  is  occurring.  A  narrowband  signal  as  well  as  a  narrowband 
jammer  will  be  used  to  facilitate  the  analysis  of  beamformer  behavior. 

A  particularly  simple  beamformer  structure  and  signal/jammer  scenario  were 
sought  to  provide  a  test  case  for  the  study  of  cancellation  effects  in  the  Frost 
ABF.  The  selected  problem  is  illustrated  in  Figure  6,  which  shows  a  sinusoidal 
desired  signal  S  (frequency  fs  =  0.25)  incident  on  a  two-element  array  from 
broadside  and  a  sinusoidal  jammer  J  (frequency  fj  =  0.26)  incident  from  45°. 
Element  spacing  was  a  half  wavelength  at  the  signal  frequency.  The  antenna 
elements  were  tied  to  a  four-weight  Frost  ABF  that  was  steered  to  broadside  and 
constrained  to  provide  unit  gain  and  flat  frequency  response  in  the  look  direction. 
(The  constraints  were  f  \  =  1  and  / 2  =  0,  as  indicated  in  the  equivalent  look- 
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direction  filter.)  The  weights  were  updated  using  the  constrained  LMS  algorithm 
given  in  Equations  (2-19). 

Just  enough  degrees  of  freedom  were  provided  in  this  test  case  to  allow  the 
objectives  of  signal  reception  and  jammer  rejection  to  be  realized.  Constant  gain 
in  the  signal  direction  was  assured  by  the  hard  constraints  of  the  Frost 
algorithm: 


!w ,  +  w2  =  l 

w3  +  w4  =  0 


(2-20) 


The  constraints  consumed  two  of  the  four  available  degrees  of  freedom;  the 
remaining  two  degrees  of  freedom  were  adequate  for  eliminating  the  jammer. 

The  requirements  for  eliminating  jammer  energy  can  be  readily  derived  for 
the  simple  case  at  hand.  The  array  output  attributable  to  the  jammer  is  given 
by 

yj(t)=  «;u"‘(wi-/w3]  +  e;(u',‘^")(w2-j;'w4]  ,  (2-21) 


where  <j>et  is  the  element-to-element  phase  shift  for  the  jammer  signal.  The 
element- to- element  phase  shift  may  be  computed  from  the  array  geometry  and 
the  signal  and  jammer  parameters: 


2k 


=  2.31  radians 


(2-22) 


Inserting  the  value  for  <fret,  expanding,  and  equating  real  and  imaginary  parts  to 


zero  yields 


{w |  -  .6739  w 2  -  .7388  w4  =  0 

-  .7388  w2 - w3  +  .6739  w4  =  0  *  (2‘23) 

The  equations  2-20  arising  from  the  constraints  and  the  equations  just  given 
may  be  solved  simultaneously  to  arrive  at  the  ideal  weight  values  for  signal 
reception  and  jammer  nulling.  The  desired  weights  are: 


Wj  =  0.50 
v/ 2  =  0.50 
w3  =  -0.22 
w4  =  0.22 


(2-24) 


These  are  the  Wiener  weight  values  that  will  be  approached  as  the 
adaptation  constant  n  is  decreased  to  small  values  and  adaptation  proceeds  very 
slowly.  Ideal  beamformer  operation  (i.e.,  perfect  signal  reception  and  jammer 
nulling)  is  obtained  with  these  weights.  In  practice,  however,  dynamic  signal 
environments  may  make  high  adaptation  speeds  necessary,  and  weight  values 
may  be  far  from  the  Wiener  ideal.  It  is  of  interest  to  examine  weight  behavior 
under  these  more  demanding  circumstances. 


Figure  7  shows  the  weight  dynamics  over  a  period  of  256  adaptation  cycles 
for  the  test  case  outlined  above  and  /i  =  0.1.  The  epoch  shown  is  after  initial 
transients  have  died  away;  the  non-Wiener  behavior  exhibited  here  will  continue 
indefinitely.  Each  weight  can  be  seen  to  oscillate  about  the  desired  Wiener 
weights;  i.e.,  the  average  values  of  the  weights  are  the  Wiener  weights.  The 
oscillation  frequency  of  0.01  /,im p  is  the  difference  frequency  between  the  desired 


Figure  7.  Weight  Dynamics  in  the  Frost  ABF. 


This  is  the  same 


signal  at  0.25  f,amp  and  the  jammer  signal  at  0.26  f,amp  . 
weight  behavior  noted  by  Glover  [11,12]  in  a  related  noise-cancelling  problem. 

Figure  8  documents  the  behavior  of  the  weights  when  the  simulation 
experiment  is  altered  by  abruptly  switching  off  the  signal  at  the  time  indicated. 
The  oscillatory  behavior  of  the  weights  ceases,  and,  after  a  brief  transient  period, 
the  weights  reach  and  hold  the  optimal  values  that  were  calculated  earlier.  It 
thus  appears  that  signal  energy  is  needed  to  support  the  non-Wiener  behavior 
that  has  been  observed.  The  precise  role  of  the  signal  in  the  adaptive  system 
may  be  better  understood  by  examining  the  interactions  of  the  desired  and 
jammer  signals  and  the  dynamic  weights. 

From  Figure  7  it  is  clear  that  during  signal  presence  each  of  the  weights  is  of 
the  form 

y»i  =  C{  +  Af  «n(wAf  +  <t>i)  ,  (2-25) 

where  C{  is  a  constant  giving  the  mean  weight  value,  cja  =  (2tt)(0.01  f,amp)  is 
the  radian  difference  frequency,  and  A,-  and  <f>{  are  constants  expressing  the 
amplitude  and  phase  of  the  sinusoidal  component  of  the  weight.  Arbitrarily 
selecting  w  i  as  the  phase  reference  and  measuring  the  remaining  constants  yields 
the  following  expressions  for  the  weights: 
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Wj  =  0.50+  0.46  sin 

w2  =  0.50  +  0.46  sin  {u&t-<f>0-n) 

w3  =  -0.22  +  0.46  sin  (wA«-^-ff/2)  *  (2‘26) 

w4  =  0.22  +  0.46  sin  (u^t-<f>0-Sir/2) 

i 

Note  that  the  oscillations  of  weights  w  j  and  w  2  are  180°  out  of  phase  so  that  the 
hard  constraint  for  these  two  weights  is  always  satisfied.  Similarly,  weights  w3 
and  w4  are  antiphased  and  satisfy  the  constraint  at  all  times. 

For  purposes  of  analysis  the  beamformer  of  Figure  6  may  be  redrawn  to 
reflect  the  observed  nature  of  the  weight  behavior.  Figure  0  shows  a  revised 
structure  in  which  each  weight  has  been  replaced  by  the  parallel  combination  of  a 
fixed  gain  and  a  mixer.  The  fixed  gains  are  set  to  the  Wiener  weights  as  required 
by  the  observed  average  values  of  the  weights.  The  mixers  have  introduced  into 
them  a  set  of  sinusoidal  voltages  with  radian  frequency  equal  to  as  required 
by  the  observations  regarding  the  oscillatory  component  of  the  weights.  The 
relative  phases  of  the  mixer  input  voltages  are  shown  in  the  phasor  diagram  inset 
into  the  upper  right  of  the  figure. 

It  is  evident  from  Figure  9  that  the  oscillatory  behavior  of  the  weights 
creates  a  considerably  more  complicated  signal  set  in  the  adaptive  beamformer 
than  would  be  the  case  with  constant  weights.  The  fixed  gain  associated  with 
each  weight  contributes  one  component  at  the  signal  frequency  and  a  second 
component  at  the  jammer  frequency.  The  mixer,  however,  yields  sum  and 
difference  components  for  both  the  signal  and  the  jammer.  The  net  contribution 
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from  each  weight  is  six  components  at  four  different  frequencies;  the  signal  set 
entering  the  final  summing  junction  in  Figure  9  consists  of  24  different 
components.  A  detailed  examination  of  these  components  is  necessary  to 
understand  signal  cancellation  in  even  this  simple  case. 

Fortunately,  the  apparent  complexity  of  the  24-signal  set  is  greater  than  its 
actual  complexity.  First,  it  should  be  recalled  that  the  Wiener  weights  are 
known  in  this  case  to  provide  perfect  signal  reception  and  perfect  jammer  nulling. 
It  follows  that  the  eight  components  attributable  to  the  fixed  gains  combine  to 
simply  yield  a  single  component  equal  to  the  desired  signal  as  seen  at  either  array 
element.  Second,  it  should  be  noted  that  in  both  beamformer  columns  (i.e.  at 
weights  w  1  and  w2  and  at  weights  w3  and  w4)  the  desired  signal  is  cophased 
while  the  weights  are  antiphased.  All  eight  terms  arising  from  mixer  action  on 
the  desired  signal  therefore  sum  to  zero  in  the  output  summing  junction.  It 
should  be  noted  that  components  generated  from  the  desired  signal  by  weight 
dynamics  will  always  cancel  in  a  properly  steered  Frost  ABF  due  to  the  phase 
relationships  enforced  by  array  steering  and  the  hard  constraints. 

At  this  point  the  only  components  that  have  not  been  examined  are  those 
generated  by  mixer  action  on  the  jammer  signals.  The  output  signal  may  be 


written  as 


+  [cos  ujt)\A  sin (wA<  -  ^0)] 

+  [cos(<^y#  -  <t>")]  [A  sin(cjAf  -<f>0-  >r)j 
+  [sin  ujt\  [A  sin(wAf  -  -  ^r/2)] 

+  [sin(wy*  -  [A  sin(wAf  -<f>0-  3jt/2)|  ,  (2-27) 

where  A  is  the  amplitude  of  the  weight  oscillations.  Expanding  the  products  of 
signals  and  weights  yields 

y(f)  =  cos  ust 

+  y  sin[(wA  +  uj)t  -  ^ol  +  “  sin[(-w5)<  -  ^0] 

+  y  sin[(a;A  +  <jj)t  -  <f>0  -  <f>tt  -  it]  +  y  sin[(-ws)f  -  <f>0  +  <f>„  -  Jr] 

-  J  sin[(o;A  +  <jj)t  +  y  sin[(-a/5)f  -  <f>0] 

A  a 

-  y  sin[(wA  +  uj)t  -  <(>0  -  <j>te  -  n]  +  y  sin[(-w5)/  -  <j>0  +  <f>ee  -  7r] 

(2-28) 

Further  simplifications  are  obvious  at  this  point.  Cancellation  eliminates  all 
the  terms  at  the  radian  frequency  (wA  +  uij).  Additionally,  terms  at  -ujs  can  be 
combined  and  then  reexpressed  using  the  relationship 

sin(-0)  =  -  sin(0)  .  (2-29) 

These  steps  plus  inserting  the  measured  value  of  A  yield 

y(<)  =  1.0  cos  ust  -  0.46  sin(u>s<+  <f>t)  -  0.46  «m(a/5f+  <t>0-<t>ee+  ir)  .  (2-30) 

The  above  expression  for  y(t)  indicates  the  interesting  nature  of  the 
beamformer  output  in  the  case  at  hand.  Weight  dynamics  have  generated  a 


number  of  components  at  various  frequencies,  but  weight  phases  are  such  that  all 
components  except  those  at  the  signal  frequency  sum  to  zero  within  the 
beamformer.  The  net  result  of  the  weight  dynamics  is  a  pair  of  synthesized 
components  at  the  signal  frequency  that,  when  added  to  the  actual  signal,  serve 
to  drive  the  beamformer  output  power  below  the  desired  level. 

Figure  10  is  a  phasor  diagram  that  shows  the  three  output  signal 
components  and  their  sum.  The  output  signal  amplitude  is  0.64,  not  1.00  as 
desired;  output  power  is  approximately  3.9  dB  below  the  desired  level.  The 
signal  has  been  partially  cancelled  by  non-Wiener  behavior  in  the  weights  even 
though  the  Frost  constraint  is  perfectly  sustained  in  the  look  direction.  Faster 
adaptation  would  cause  even  more  signal  loss. 

E.  An  Equivalent  Adaptive-Noise-Cancelling  Problem 

The  problem  set  forth  in  Figure  6  can  be  treated  as  an  adaptive-noise- 
cancelling  problem  if  it  is  imagined  that  the  desired  and  jammer  signals  are 
available  separately.  The  equivalent  problem  structure  is  illustrated  in  Figure  11. 
The  desired  signal  may  be  taken  straight  through  to  the  summing  junction  since 
the  hard  constraints  specify  unity  gain  and  a  flat  frequency  response.  The 
jamming  signal  is  brought  to  the  summing  junction  through  an  adaptive  filter 
that  emulates  the  jammer-direction  gain  of  the  Frost  beamformer.  An  in-phase 
gain  (7/  and  a  quadrature  gain  Gq  are  able  in  this  special  case  to  give  complete 
control  over  both  the  amplitude  and  phase  of  the  jammer  contribution. 
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Expressions  relating  Gj  and  Gq  to  the  ABF  weights  wlf  w2,  w3)  and  w4 
are  needed  to  complete  the  specification  of  the  equivalent  problem.  The  required 
expressions  may  be  derived  from  the  equation  for  the  jammer  output  from  the 
ABF: 

Hj{t)  =  w1cosu>/<  +  w2cos(w/f-0ee)  +  w3sinw/<+  w4  sinfcjyt-^,,) 

(2-31) 

where  <j>ee  is  the  element-to-element  phase  delay  for  the  jammer  signal.  After 
making  use  of  trigonometric  identities  and  regrouping  terms,  the  equation 
becomes 

y/(0  =  Iwi  +  (cos  <t>tt )w 2-(sin  ^JwJ  cos  Ujt 

+  ((sin  <f>et )w 2  -f-  w3  +  (cos  ^JwJ  sin  ujt  .  (2-32) 

Comparison  of  this  equation  with  the  structure  of  Figure  11  yields  the 
relationships 

Gl  =  -(wi+  (cos  <f>te  )w  2  -(sin  ^ejwj  (2-33) 

Gq  =  -[(sin^„)w2+  w3+  (cos^jwj  . 

The  equivalent  system  involving  G\  and  Gq  allows  a  considerably  simpler 
characterization  of  the  weight  dynamics  associated  with  signal  cancellation  than 
would  have  been  possible  in  the  four-dimensional  weight  space  of  the  original 
problem.  In  particular,  the  reduced  dimensionality  simplifies  comparisons  of  the 
dynamic  solution  obtained  during  cancellation  with  the  Wiener  solution  that  is 
approached  at  low  adaptation  rates. 


The  Wiener  solution  in  terms  of  Gt  and  Gq  may  be  determined  by 


inspection  for  the  ANC  structure: 
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Thus  it  is  seen  that,  in  the  two-space  defined  by  plotting  Gj  along  the  x  axis  and 
Gq  along  the  y  axis,  the  origin  represents  the  solution  that  gives  perfect  signal 
recovery  and  perfect  jammer  rejection. 

Figure  12  is  a  plot  that  shows  the  loci  of  several  dynamic  solutions  in  the 
G 'i  -  Gq  plane.  Each  solution  is  obtained  by  repeating  the  narrowband 
simulation  of  the  previous  section  for  the  indicated  value  of  ft  and  computing  Gj 
and  Gq  from  the  weights  {w,}.  At  /i=. 001,  the  solution  locus  is  a  tiny  circle 
that  blurs  into  a  dot  at  the  origin,  i.e.,  the  solution  almost  exactly  duplicates  the 
Wiener  solution.  Signal  cancellation  is  negligible  in  this  case.  As  /i  is  increased, 
the  solution  loci  expand  to  form  circles  *  about  the  origin,  and  the  amount  of 
signal  cancellation  is  increased.  Cancellation  exceeds  90  percent  at  the  point 
where  stability  considerations  forbid  further  increases  in  /i. 

It  is  worthwhile  to  note  the  range  of  n  over  which  departure  from  the 
Wiener  solution  is  visible  in  Figure  12.  At  /i  =  0.001  the  deviation  from  the 
Wiener  solution  is  barely  perceptible.  Instability  does  not  occur  until  n  reaches 

*Some  departs!*  from  •  circa  Ur  local  ii  evident  for  (1  =  0.1.  The  effect  ii  dne  to  tecosd-order  weight  oscilla¬ 
tion  at  the  lain  freqaencjr,  i.e.,  at  f  s  +  f  /■  Glover  reported  similar  effect*  is  111). 
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beyond  1.00  (approximately  1.15).  This  gives  a  useful  range  of  ft  of  about  three 
orders  of  magnitude  over  which  departure  from  the  Wiener  solution  is  noticeable. 
This  is  a  rather  broad  operating  region  that  includes  the  values  of  ft  for  which 
convergence  is  most  rapid.  In  many  cases  it  will  be  very  unattractive  to  sacrifice 
orders  of  magnitude  in  convergence  rate  in  order  to  obtain  a  close  approximation 
to  the  Wiener  solution.  This  is,  unfortunately,  the  nature  of  the  tradeoff  that 
must  be  addressed  when  conventional  adaptive  beamformers  are  applied. 

The  loci  of  Figure  12  have  several  interesting  interpretations.  Considered 
simply  as  the  path  traced  by  the  tip  of  a  gain  vector,  each  locus  demonstrates  the 
same  cyclic  weight  behavior  that  was  shown  earlier  in  Figure  7.  The  number  of 
samples  required  to  complete  one  circuit  of  a  locus  is  100;  this  is  the  period  of  the 
difference  frequency  between  the  signal  frequency  of  0.25  f8amp  and  the  jammer 
frequency  of  0.26  ftamp . 

The  radius  of  each  locus  may  be  interpreted  as  a  measure  of  the  peak  short- 
term  correlation  between  the  desired  signal  and  the  jammer.  The  time  window 
over  which  the  short-term  correlation  is  computed  is  governed  by  the  parameter 
ft:  large  values  of  ft  lead  to  rapid  adaptation  and,  consequently,  narrow  time 
windows.  Peak  correlation  between  the  desired  signal  and  the  jamming  signal 
can  be  quite  high  if  only  a  narrow  time  window  is  considered.  The  peak  short¬ 
term  correlation  for  ft  =  1.0,  for  example,  is  almost  unity.  The  weight  dynamics 
serve  to  create  a  modified  jamming  signal  in  which  this  short-term  correlation  is 
maintained  near  its  peak  value.  Appendix  A  considers  the  dynamic  solution  from 
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the  viewpoint  of  short-term  correlations  and  develops  a  closed-form  description  of 
weight  behavior  when  both  signal  and  jammer  are  sinusoidal. 

Another  way  to  view  the  loci  of  Figure  12  is  to  treat  them  as  traces  of  the 
minima  of  dynamic  error  surfaces.  It  is  well  known  that,  in  the  stationary  case, 
the  Wiener  solution  defines  the  minimum  of  a  quadratic  error  surface.  In  the 
case  at  hand,  the  error  surface  may  be  treated  as  static  only  at  small  values  of  ft. 
The  minimum  of  the  surface  in  the  stationary  case  is  not  at  zero  error,  but  is 
elevated  above  zero  by  the  power  level  of  the  desired  signal.  At  low  values  of  ft 
the  Wiener  solution  is  closely  approximated,  and  the  desired  signal  is  delivered 
essentially  without  loss  at  the  beamformer  output.  As  p  is  increased,  the  error 
surface  must  be  treated  as  dynamic.  The  minimum  of  the  error  surface  no  longer 
coincides  with  the  Wiener  solution.  Instead,  for  the  simple  case  under 
consideration,  the  minimum  lies  on  a  circle  centered  on  the  Wiener  solution. 
Furthermore,  the  dynamic  minimum  does  not  lie  at  the  same  error  level  as  the 
Wiener  solution.  The  elevation  of  the  minimum  is  decreased  by  the  amount  of 
signal  cancellation  attainable  at  a  given  value  of  /*.  The  decrease,  expressed  as  a 
fraction  of  the  desired  signal  power,  is  related  to  the  quantity  termed 
miaadjuatmcnt  [3] 

misadjustment  =  A/(p)  —  ^ll/  (01 — ^f(y  (Offl.  (2-35) 

where  y*{t)  is  the  output  with  the  optimal  weight  vector  W*.  Because  £[y2(/)J 
is  less  than  £((?"(*  ))*l»  misadjustment  is  a  negative  quantity  in  the  region  of 
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interest.  It  is,  however,  the  magnitude  of  misadjustment  that  is  important,  not 
the  sign.  The  contours  of  Figure  12  have  been  labelled  to  show  \M(p)  |  levels  as 
well  as  p  values. 

The  discussion  of  the  error  surface  that  has  just  been  given  differs  somewhat 
from  earlier  treatments.  For  a  given  ABF  structure,  the  error  surface  has 
previously  been  considered  simply  as  a  function  of  the  input  signals;  changes  in 
the  error  surface  were  attributable  to  changes  in  the  input  signal  description. 
Here  the  error  surface  has  been  described  as  a  function  of  the  adaptation 
parameter.  In  other  words,  the  error  surface  must  be  described  in  terms  of  the 
signal  environment  perceived  by  the  beamformer,  not  simply  in  terms  of  some 
detached  statistical  characterization  of  the  incident  signals. 

The  view  of  the  solution  set  that  has  just  been  given  indicates  the  nature  of 
the  problem  at  hand.  There  is  a  direct  coupling  between  the  adaptation 
parameter  that  is  chosen  and  the  solution  that  is  delivered  by  the  ABF.  Low 
values  of  the  adaptation  parameter  p  yield  a  solution  that  conforms  closely  to  the 
Wiener  solution  and  satisfies  the  requirement  for  signal  preservation.  These 
values  of  p  fail,  however,  to  meet  the  requirements  of  responding  rapidly  and 
tracking  a  dynamic  signal  environment.  Higher  values  of  p  provide  better 
tracking  of  changes  in  the  environment,  but  the  non-Wiener  solutions  involve 
sacrifices  in  the  quality  of  the  recovered  signal.  Given  the  beamformer  under 
consideration,  nothing  can  be  done  beyond  striking  the  best  compromise  between 
signal  quality  and  tracking  capability.  A  change  in  the  nature  of  the  beamformer 
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m.  A  COMPOSITE  BEAMFORMER  THAT  ELIMINATES  SIGNAL 
CANCELLATION 

Chapter  2  has  illustrated  the  signal  cancellation  effects  that  can  arise  in 
relatively  simple  signal/jammer  scenarios  and  has  demonstrated  that  even  the 
most  rigid  constraints  can  fail  to  preserve  the  desired  signal.  This  chapter  turns 
from  a  discussion  of  the  problem  to  the  description  of  a  solution  and  shows  how 
signal  cancellation  can  be  avoided  at  the  price  of  some  increase  in  beamformer 
complexity. 

Chapter  3  consists  of  three  sections.  The  first  section  builds  upon  the 
background  provided  by  Chapter  2  and  reformulates  the  adaptive  beamformer 
problem  in  such  a  way  that  jammer  nulling  can  be  accomplished  at  high 
adaptation  rates  without  signal  cancellation.  The  beamformer  structure  that 
arises  out  of  the  reformulation  is  termed  a  composite  beamformer  or  CBF.  The 
next  section  discusses  the  key  issue  of  signal  relationships  within  the  new 
beamformer.  A  demonstration  of  the  performance  improvement  afforded  by  the 
new  beamformer  is  given  in  the  third  section  by  revisiting  the  wideband  problem 
from  Chapter  2. 

A.  Problem  Reformulation 

Chapter  2  and  Appendix  A  described  the  differing  solutions  derived  by  a 
hard-constrained  beamformer  at  various  adaptation  rates.  The  solutions  differ 
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because,  from  the  beamformer  viewpoint,  the  problem  changes.  Unfortunately, 
the  viewpoints  of  the  beamformer  and  the  system  designer  begin  to  clash  as 
signal  cancellation  becomes  significant.  The  designer  is  not  interested  in  the 
minimization  of  output  power  at  all  costs,  and  will  typically  invoke  constraints  in 
an  attempt  to  restrict  the  minimization  process.  The  beamformer,  however,  is 
designed  to  relentlessly  pursue  minimization.  At  high  adaptation  rates  the 
beamformer  weights  possess  the  mobility  to  exploit  short-term  correlations 
between  signal  and  jammers  and  thereby  circumvent  the  constraints.  Design 
objectives  are  simultaneously  bypassed,  and  system  performance  is 
unsatisfactory.  A  problem  reformulation  is  needed  that  harmonizes  design 
objectives  with  the  realities  of  beamformer  behavior  at  high  adaptation  rates;  this 
section  pursues  that  reformulation. 

Two  observations  can  be  made  at  this  point  that  are  useful  in  developing  a 
problem  reformulation.  One  observation  b  that  interaction  between  the  desired 
signal  and  the  jammer  b  the  root  of  the  cancellation  phenomenon.  The  nature  of 
thb  interaction  was  demonstrated  in  some  detail  in  Chapter  2  for  the  case 
involving  a  signal  and  a  jammer  that  are  narrowband.  It  was  shown 
experimentally  and  analytically  that  the  presence  of  both  signal  and  jammer 
energy  b  a  prerequbite  for  signal  cancellation.  In  particular,  it  was  shown  that 
the  output  signal  b  the  “target”  of  the  cancellation  process  and  that  short-term 
correlation  between  signal  and  jammer  waveforms  b  the  phenomenon  that  makes 
cancellation  possible. 
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A  second  observation  that  can  be  made  is  that  the  signal  plays  no  role  in  the 
Wiener-solution  calculation  in  a  perfectly  steered  Frost  ABF.  In  other  words,  it 
is  completely  equivalent  to  write 

W  =  R^C\CT  Rj„  C ]  L  (3-1) 


rather  than  the  earlier  statement 


W  =  Rj&C\CT  Rj&C]l  , 


(2-13) 


where  Rxx  =  Rss  +  RNff  as  originally  assumed.  This  point  can  be  appreciated 
intuitively  from  the  fact  that  the  look-direction  response  is  determined 
exclusively  by  the  hard  constraints,  not  by  the  desired  signal.  The  equivalence  of 
(3-1)  and  (2-13)  may  be  rigorously  demonstrated  by  substituting  for  RXx  in  (2- 
13)  and  simplifying.  An  appropriate  substitution  is 


Rxx  =  Rss  +  Rnn  —  CT  rrr  C+  R, 


where  Rrr  is  the  autocorrelation  matrix  for  the  vector  of  signal  voltages 
appearing  at  the  taps  of  any  one  of  the  beamformer  filters.  Appendix  B  traces 
the  somewhat  lengthy  proof  of  the  equivalence  for  both  a  narrowband  processor 
using  a  single  complex  weight  per  beamformer  channel  and  a  wideband  processor 
using  tapped- delay-line  filters  for  each  beamformer  channel.  Yet  another  way  to 
explore  the  role  of  the  signal  is  to  consider  the  problem  in  the  context  of  an 
equivalent  generalized  sidelobe  canceller;  Appendix  C  discusses  this  approach  to 
the  problem. 
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Reconsideration  of  the  adaptive  beamformer  problem  has  thus  far  indicated 
1)  that  the  signal  drives  the  process  of  signal  cancellation  that  results  in  its  own 
demise  and  2)  that  the  signal  has  no  role  in  the  determination  of  a  set  of  weights 
that  optimizes  array  performance.  These  two  points  argue  for  exclusion  of  the 
desired  signal  from  the  beamformer.  The  overlooked  consideration  is,  of  course, 
that  the  original  objective  was  to  recover  the  signal.  The  remaining  step  is  to 
harmonize  the  objectives  of  cancellation-free  adaptation  and  successful  signal 
reception  by  devising  a  structure  that  excludes  the  signal  from  the  adaptive 
process  but  allows  the  signal  to  pass  through  to  the  system  output.  The  problem 
reformulation  thus  amounts  to  more  precisely  defining  the  role  of  the  desired 
signal  in  an  adaptive  array. 

Figure  13  illustrates  a  beamformer  based  upon  the  problem  reformulation. 
The  array  has  been  augmented  so  that  subarrays  consisting  of  multiple  elements 
appear  in  place  of  the  individual  elements  of  the  original  array.  A  preprocessor 
operates  upon  the  received  signals  from  the  augmented  array  to  generate  an 
environment  that  is  free  of  desired-signal  content.  A  Frost  ABF  operates  in  this 
synthetic  environment  and  derives  weights  that  are  copied  to  a  slaved 
beamformer.  The  slaved  beamformer  has  the  same  signal- path  structure  as  the 
adaptive  beamformer,  but  is  connected  directly  to  selected  antenna  elements.* 
Desired-signal  components  are  present  in  the  slaved  beamformer,  and  a  useful 

*lt  it  alio  possible  to  connect  the  staved  beamformer  to  the  snbarrayi  through  a  preprocessor  that  makes  use  of 
meltiple  elements  of  each  snbarray.  The  pattern  response  of  the  element  combination  generated  by  this  second  preproces¬ 
sor  most  have  nonsero  response  in  the  look  direction. 
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output  may  be  drawn  from  it. 


Taken  as  a  whole,  the  structure  shown  in  Figure  13  is  termed  a  composite 
beamformer  (CBF).  The  key  elements  are  1)  an  augmented  array,  2)  a 
preprocessor  that  excludes  the  desired  signal  from  the  adaptive  process,  3)  an 
adaptive  beamformer  that  can  be  constrained  to  control  look-direction  response 
while  nulling  jamming  signals,  and  4)  a  slaved  beamformer  that  is  used  to 
implement  the  computed  solution  and  recover  the  desired  signal.  Three  of  these 
elements,  the  array,  the  preprocessor  and  the  adaptive  beamformer,  afford 
considerable  flexibility  in  that  a  variety  of  specific  realizations  are  possible;  the 
slaved  beamformer  design  is  inflexible  in  the  sense  that  it  mirrors  the  adaptive 
beamformer  design.  Appendix  D  describes  an  alternative  CBF  realization  based 
upon  Widrow’s  pilot-signal  algorithm;  other  realizations  may  also  be  derived. 

B.  Signal  Relationships  In  the  CBF  Master  and  Slave 

Because  CBF  operation  involves  the  derivation  of  weights  in  the  master 
beamformer  for  application  in  the  slave,  it  is  essential  that  the  signal  sets  in  the 
two  beamformers  be  closely  related.  Ideally,  the  two  signal  sets  would  be 
identical,  save  for  the  absence  of  look-direction  components  in  the  master.  This 
would  assure  that  weights  derived  in  the  master  were  equally  appropriate  for  and 
effective  in  the  slave.  Unfortunately,  it  is  not  possible,  in  general,  to  attain  the 
ideal  signal  relationships  between  the  two  beamformers,  and  compromises  must 
be  allowed.  This  section  begins  a  discussion  of  the  nature  of  the  compromises 
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that  are  necessary  and  of  the  consequences  of  those  compromises. 

The  basic  difficulty  in  attaining  the  ideal  signal  sets  lies  in  removing  the 
look-direction  signal  from  the  master  beamformer  without  perturbing  the  various 
noise  signals.  If  an  omnidirectional  element  is  used  for  the  slaved  beamformer, 
then  the  perfect  subarray /preprocessor  combination  serving  the  master  must 
have  the  pattern  response  shown  in  Figure  14(a).  The  omnidirectional  pattern  is 
perfectly  replicated,  except  for  an  infinitely  deep,  vanishingly  narrow  notch  in 
precisely  the  look  direction.  Obviously,  a  pattern  of  this  nature  can  neither  be 
realized  nor  successfully  applied;  some  modification  of  interference  signals 
arriving  near  the  look  direction  must  be  accepted  if  the  desired  signal  is  to  be 
excluded  based  upon  directional  information. 

A  simple  means  of  realizing  a  look-direction  null  is  to  use  a  two-element 
subarray,  apply  any  necessary  steering  delays  to  cophase  the  desired  signal,  and 
difference  the  element  signals.  Figure  14(b)  gives  the  response  of  such  a  subarray 
at  the  frequency  for  which  the  elements  are  separated  by  one-half  wavelength. 
The  null  is  steered  to  broadside  in  this  instance. 

The  response  of  the  two-element  subarray  departs  significantly  from  the 
ideal  response.  At  ±  90®,  the  two-element  subarray  exhibits  6  dB  of  gain  over 
isotopic;  0  dB  is  at  about  ±  20®;  gains  within  ±  10®  of  the  look  direction  fall 
seriously  below  0  dB.  The  two-element  pattern  does,  however,  exhibit  the  key 
features  of  a  deep,  wideband  null  in  the  look  direction  and  reasonable  gains  at 
angles  removed  from  the  look  direction.  These  features  plus  the  great  simplicity 
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of  the  subarray /preprocessor  make  the  two-element  subarray  worthy  of  more 
detailed  consideration. 

The  general  CBF  structure  of  Figure  13  may  be  specialized  to  yield  the  CBF 
shown  in  Figure  15.  Here  the  desired  signal  is  assumed  to  be  incident  from 
broadside  so  beam  steering  may  be  neglected.  The  preprocessor  is  realized  by 
using  two-element  subarrays  in  a  simple  element-differencing  scheme.  The 
sharing  of  elements  between  subarrays  provides  very  efficient  use  of  elements; 
only  a  single  element  is  needed  beyond  those  ordinarily  required  for  a  comparable 
array.* 

A  Frost  ABF  operates  upon  the  preprocessor  output  to  derive  a  set  of 
weights  that  minimizes  error  power  subject  to  a  set  of  look-direct  constraints.  It 
should  be  noted  that,  in  contrast  to  conventional  adaptive  beamformers,  the  error 
power  is  distinct  from  the  CBF  output  signal  power.  Error  power  can  be  driven 
to  zero  without  endangering  the  desired  signal  component  in  the  output.  This  is, 
in  fact,  the  ideal  state  of  affairs. 

Because  weights  derived  in  the  Frost  ABF  will  be  applied  in  the  slaved 
beamformer,  it  is  vital  that  relative  signal  phases  match  in  the  two  beamformers. 
The  uniform  structure  that  exists  in  the  preprocessor  assures  that  this  phase 
matching  will  be  obtained.  The  set  of  two-element  subarrays  feeding  the  Frost 
ABF  echoes  the  structure  of  the  array  of  omnidirectional  elements  serving  the 

*Elem*Dt  ibaring  among  enbarrays  does  introdoce  thermal-noiie  corrclatiom  is  the  master  beamformer  that  are 
aot  present  in  the  tlare.  Generally,  however,  thermal  noiee  powers  will  be  relatively  low,  and  thee#  epnrione  correlations 
will  not  eerionely  affect  performance. 


Jammer 


Figure  15.  Frost-based  CBF  with  a  Subtractive  Preprocessor  for  Signal  Nulling 


slaved  beamformer. 


Phase  relationships  between  signals  in  the  ABF  and  the  slaved  beamformers 
are  made  clearer  in  Figure  16.  The  jammer  components  received  by  the 
omnidirectional  elements  are  indicated  by  a  set  of  equal-amplitude,  uniformly 
spaced  phasors  J j,  J2,  ^z>  •••»  ^K-i*  *^K-  The  preprocessor  operates  upon  these 
inputs  to  produce  the  phasor  outputs  •••>  and  These 

phasors  are  also  equal  in  amplitude  and  have  the  same  phase-angle  separations  as 
the  received  jammer  components.  The  preservation  of  relative  phase  in  the 
preprocessor  assures  that  weights  that  generate  a  null  in  the  ABF  will  also 
generate  a  null  in  the  slaved  beamformer. 

The  phasor  argument  as  advanced  in  Figure  16  applies  to  a  single  jammer  at 
a  single  frequency.  Linearity  and  superposition  apply,  however,  and  show  that 
phase  relationships  are  preserved  for  multiple  jammers  and  for  broadband  as  well 
as  narrowband  signals. 

The  uniform  linear  array  provides  an  attractive  structure  for  the  CBF 
because  there  is  the  option  of  element  sharing  between  subarrays.  A  regular 
array  structure  is  not,  however,  a  prerequisite  for  the  CBF.  The  fundamental 
requirement  is  for  phase  matching  between  the  master  and  slave  beamformers. 
Phase  matching  may  be  obtained  for  an  arbitrary  array  geometry  by  augmenting 
each  original  element  to  form  identical  subarrays  at  the  element  locations. 
Identical  preprocessors  may  then  be  used  to  form  subarray  responses  with  nulls  in 
the  selected  look  direction.  The  preprocessor  outputs  must  be  cophased  for  the 


look  direction  and  applied  to  the  adaptive  beamformer  while  the  original  element 
outputs  are  cophased  and  applied  to  the  slave. 

C.  CBF  Performance  —  An  Example 

The  CBF  described  in  the  previous  section  is  capable  of  delivering  much 
better  performance  at  high  adaptation  speeds  than,  say,  a  comparable  Frost 
ABF.  The  performance  contrast  may  be  conveniently  illustrated  by  returning  to 
the  wideband  problem  of  Chapter  2  for  which  the  Frost  ABF  performance  has 
already  been  demonstrated. 

Figure  17  diagrams  the  wideband  problem  with  a  Frost-based  CBF  replacing 
the  original  Frost  ABF.  The  signal  and  jammer  descriptions  are  identical,  as  are 
the  look-direction  constraints.  A  third  antenna  element  has  been  added  in  order 
to  provide  appropriate  inputs  for  the  subtractive  preprocessor.  Two  of  the 
elements  also  serve  the  slaved  beamformer,  which  makes  use  of  weights  from  and 
functions  in  parallel  with  the  adaptive  beamformer. 

The  antenna  pattern  at  the  jammer  frequency  of  0.25  /Mmp  is  shown  in  Part 
a)  of  Figure  18.  This  plot  was  generated  after  initial  transients  had  died  away 
and  jammer  nulling  was  essentially  complete.  A  deep  notch  exists  at  the  jammer 
arrival  angle  of  30*. 

Parts  b)  and  c)  of  Figure  18  show  the  frequency  response  in  the  jammer 
direction  and  look  direction.  The  jammer-direction  plot  reveals  a  deep  notch  at 
the  jammer  frequency.  The  look-direction  plot  confirms  that  the  unit-gain,  all- 
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pass  response  specified  by  the  constraints  has  been  attained. 

Some  of  the  contrasts  between  Frost  ABF  and  CBF  performances  may  be 
seen  by  comparing  Figure  3  and  Figure  18.  The  notch  depths  in  both  Part  a) 
and  Part  b)  of  the  figures  are  greater  for  the  CBF  than  for  the  Frost  ABF,  thus 
indicating  improved  jammer  rejection. 

A  more  dramatic  view  of  the  performance  contrasts  is  provided  by  the 
spectral  plots  of  Figure  19.  Parts  a)  and  b)  of  the  figure  show  ensemble  averages 
of  the  desired  signal  spectrum  and  the  jammer  spectrum.  Part  c)  repeats  the 
ensemble-averaged  Frost  output  spectrum  that  was  shown  in  Figure  4.  Part  d) 
provides  the  ensemble-averaged  CBF  output  spectrum.  The  signal-cancellation 
effects  seen  in  Part  c)  are  not  present  in  the  CBF  spectrum,  and  extremely  close 
matching  with  the  desired-signal  spectrum  is  evident. 

The  differences  between  Frost  ABF  and  CBF  performance  may  also  be 
clearly  seen  in  the  time  domain.  Figure  20  contains  256-sample  segments  of  time 
domain  data  taken  from  the  beginning  of  the  simulation.  Part  a)  illustrates  the 
desired-signal  waveform;  Part  b)  the  Frost- ABF  output;  and  Part  c)  the  CBF 
output.  Both  output  waveforms  show  the  transient  associated  with  jammer 
nulling.  Nulling  should  be  complete  within  50  samples,  and  good  tracking  of  the 
desired  signal  should  be  seen.  The  Frost-ABF  output,  however,  shows 
considerable  distortion  due  to  signal  cancellation.  CBF  tracking,  on  the  other 
hand,  is  essentially  perfect  after  jammer  nulling. 


Input  and  Output 

a)  Signal 

b)  Jammer 

c)  ABF  Output 

d)  CBF  Output 


JO 


0.5 


ju  =  1.0  x  10' 
0.5 


jj  «  1.0  x  10“ 
0.5 


Frequency 


Frequency 


7\ 

Frequency 


for  the  Frost  ABF  and  the  Frost-based  CBF 


Figure  20.  Comparison  of  ABF  and  CBF  Output  Signal 


IV.  PERFORMANCE  OF  THE  COMPOSITE  BEAMFORMER 


The  preceding  chapter  introduced  the  composite  beamformer  (CBF)  and 
illustrated  the  performance  of  a  particular  CBF  realization  based  upon  a  simple 
subtractive  preprocessor  and  Frost’s  constrained  LMS  algorithm.  An  argument 
was  advanced  that,  due  to  the  deliberate  phase  matching  between  the  master  and 
slave  beamformers,  jammer  nulling  is  quite  good  in  the  CBF.  Additionally,  signal 
cancellation  of  the  type  that  can  be  so  destructive  in  a  hard-constrained  ABF  is 
effectively  suppressed.  A  sample  case  was  shown  in  which  the  expected 
improvements  over  a  Frost  ABF  were,  in  fact,  realized  by  a  CBF  using  a  simple 
signal-nulling  scheme. 

This  chapter  continues  the  exploration  of  CBF  performance.  The  first 
section  examines  the  issue  of  CBF  optimality  in  greater  depth.  The  second 
section  then  describes  further  simulation  experiments  designed  explicitly  to  probe 
CBF  behavior  with  regard  to  optimality.  A  third  section  considers  convergence 
time  constants  in  the  CBF  and  draws  comparisons  with  the  Frost  ABF.  The 
fourth  section  describes  CBF  weight  behavior  near  convergence  and  shows 
contrasts  between  CBF  and  ABF  behavior  in  that  important  region. 

A.  CBF  Optimality  in  the  Narrowband  Case 

The  experiment  of  the  previous  chapter  has  indicated  that,  at  least  in  a 
selected  case,  the  performance  of  the  CBF  at  high  adaptation  speeds  is  greatly 
superior  to  that  of  a  comparable  ABF.  Questions  remain,  however,  about  the 


range  of  conditions  for  which  optimal  performance  is  approached.  In  this  section 
the  issue  of  optimality  will  be  investigated  in  greater  depth  for  the  important 
case  in  which  jammer  signals  are  narrowband. 


The  discussion  in  this  section  will  center  upon  the  noise  autocorrelation 
matrix  R^n-  As  was  shown  in  Chapter  HI,  Rjw  governs  the  optimal^olution  for 
the  Frost  ABF.  Ideally,  the  adaptive  (master)  processor  in  the  CBF  would 
receive  unperturbed  noise  voltages  and  would  thus  tend  toward  the  optimal 
solution  as  indicated  by  Rtf#.  In  practice,  the  phases  of  interference  signals  are 
preserved  through  appropriate  preprocessor  design,  but  the  interference 
amplitudes  (or,  equivalently,  the  interference- to-thermal-noise  ratios)  are  altered. 
It  is  of  interest  to  consider  the  influence  of  perturbations  in  interference-to- 
thermal-noise  ratio  upon  Rfw. 

The  signal  vector  14  associated  with  the  »tA  narrowband  interference  source 
may  be  written  as 
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(4-1) 


where  /(tf,*,Ff)  is  the  pattern  response  of  the  K  identical  elements  (or  subarrays) 
constituting  the  array,  0,  is  the  arrival  angle  of  the  signal  with  respect  to  a 
reference  direction,  P{  is  the  signal  power,  J2,-  is  a  propagation  vector,  and  <f>ti  b 
the  signal  phase  at  the  /**  element  with  respect  to  a  coordinate  origin.  For  the 


important  case  of  a  uniform  linear  array  with  element  spacing  d ,  (4-1)  may  be 
specialized  to 


Mi  =  mvP.) 
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where  the  first  element  is  taken  as  the  phase  reference  and 


=  2"i~) sin  9. 


(4-3) 


is  the  element-to-element  phase  shift  for  an  arrival  angle  0{  (measured  clockwise 
from  array  broadside)  and  an  interference  wavelength  of  X,-. 

Pattern  responses  of  two  types  will  be  of  interest.  For  an  ideal 
omnidirectional  element, 

f  lie, Pi)  =  Pi  .  (4-4) 

If  two  such  elements  are  placed  a  half  wavelength  apart  and  their  outputs  are 
differenced  to  form  the  simple  nullformer  discussed  previously,  then 

fi {e,P{)  —  P,  [2  -  2  cos(?r  sin  0)]  .  (4-5) 


The  noise  autocorrelation  matrix  RN^  for  a  K-element  narrowband  array 
takes  a  simple  form  under  the  assumptions  that  the  thermal  noise  voltages  from 
the  elements  are  zero-mean  Gaussian  and  mutually  uncorrelated  and  that  the  M 
incident  narrowband  signals  have  carrier  phases  that  are  uniformly  distributed  on 


(0,27r)  and  statistically  independent  both  of  one  another  and  of  the  thermal  noises 
of  the  elements.  Specifically, 


&NN  ~  [I  +  S  Pi  $i  »  (4-6) 

»—i 

where  /  is  the  K  X  K  identity  matrix,  p,-  is  the  power  ratio  between  the  i‘* 
interference  signal  and  the  thermal  noise  voltage,  and  a2  is  the  thermal  noise 
variance.  The  power  ratio  p,-  incorporates  the  pattern  response  /(0,-,P,)  and  a2: 
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(4-7) 


Gupta  and  Ksienski  [20]  have  observed  that  correlation  matrices  having  the 
structure  of  (4-6)  may  be  inverted  into  a  convenient  form  when  the  matrices  are 
reexpressed  in  terms  of  an  orthonormal  basis  for  {14  }.  Their  method  will  be  used 
to  derive  an  expression  for  Pjyjv-  This  expression  will  then  be  used  to  determine 
conditions  under  which  the  CBF  solution  approaches  that  of  the  Frost  ABF. 

The  Gram-Schmidt  process  may  be  applied  to  construct  an  orthonormal 
basis  {$;•},  »=1,  2,  ...,  M,  for  the  {14},  *  =1,  2,  ...,  M.  As  a  preliminary  step  the 
first  basis  vector  is  computed: 


£1  = 


(4-8) 


Then  a  set  of  orthogonal  vectors  {&•},  »'=1,  2,  ...,  M,  is  determined  from 
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ft  =  ft  ~  E  U  V)  Sj 

i-i 


and  normalized  according  to 
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ft  ft 
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(4-10) 


The  vectors  {ft-}  may  now  be  written  in  terms  of  the  basis  vectors.  The 


expression  given  in  (4-0)  may  be  rearranged  to  yield 


Si  =  E  (xfSjlSj  +  ft 
l 


(4-11) 


Forming  the  transpose  in  (4-11)  and  postmultiplying  by  ft*  then  gives 


ft7**  =  a7 *’  =  (ft7*/)*7*"  =  ft7*.-*  . 


(4-12) 


where  the  orthonormality  of  {ft  }  has  been  used  to  simplify  the  expressions.  By 


using  (4-11)  and  (4-12)  together  with  the  deGnition 


0,7  4  ftV  , 


(4-13) 


ft-  may  be  written  as 


ft  =  E 

/- 1 


(4-14) 


The  expression  (4-6)  for  RNN  now  becomes 


v*.  t_*  v*  *,i  v”  ^.*i  «  ‘  -  " _ m  "w*  . . 
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(  M  [  i  1* 

E  °.*&1  \°2  (*-15) 

RNN  =  W  +  E  Pi  1  E  a*j  Uj  I 

(  i-1  u— 1  J 

L*-»  J  i 

Because  the  optimal  weights  for  both  the  Frost  ABF  and  the  Frost-based  CBF 
depend  upon  Rfi jy,  it  is  necessary  to  invert  the  expression  on  the  right-hand  side 
of  (4-15).  The  inversion  technique  depends  upon  the  matrix  inversion  lemma* 
and  is  most  readily  applied  by  first  considering  the  case  M  =  1.  In  that  case 
(4-15)  becomes 

rNN  =  V  +  Pi  aii  Sli  1  af  }  o2  (4-16) 

Application  of  the  inversion  lemma  to  (4-16)  then  yields 

Rm  ~  U  ~  Pi  Qii  4i(l  +  Sli  Pi  ah  Sli)'1  Hi  }  »  t4-17) 


which  rearranges  to  the  form 
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By  orthonormality, 


Sli  Sli  =  i 


(4-10) 


Also 
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The  femme  is  discassed  ia  Appendix  B  ia  another  context. 
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The  expression  for  in  (4-18)  therefore  becomes  simply 
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It  is  now  clear  that,  provided  px  »  1, 


=  V-alaT)^ 


(4-22) 


That  is,  if  the  jammer-to- thermal-noise  ratio  is  large,  the  inverse  of  the  noise 
autocorrelation  matrix  is  essentially  independent  of  the  precise  value  of  that 
ratio.  Because  the  optimal  solution  is  governed  by  /? jvJv>  the  solutions  toward 
which  the  ABF  and  CBF  tend  will  be  indistinguishable,  provided  the  jammer-to- 
thermal-noise  ratio  is  large  in  both  adaptive  processors. 

In  the  ABF  with  omnidirectional  elements  the  jammer-to-thermal-noise  ratio 
is  independent  of  jammer  angle,  and  the  ratio  is  thus  fixed  solely  by  relative 
power  considerations.  In  the  CBF  with  element  differencing  the  jammer-to- 
thermal-noise  ratio  is  dependent  upon  jammer  arrival  angle  in  the  manner  shown 
in  (4-5).  At  angles  approaching  0®  the  received  jammer  power  goes  toward  zero, 
and  the  transition  from  (4-21)  to  (4-22)  becomes  questionable.  It  thus  becomes 
clearer  that  ABF  and  CBF  equivalence  is  not  universal  and  that  the  nullformer 
becomes  the  crucial  factor  in  CBF  convergence  toward  the  optimal  weights. 


i^v,v,vAV.',VAv.*y,v. '-V-Y*V- . 


If  consideration  is  now  given  to  the  two-jammer  case  (4-15)  becomes 


RnN  —  {^  +  PlallAlMi  +  P2a2i&l(lT  +  P2a2lQ22SlH2 

+  P2a22a2l4a8 2  +  P2a22&2Jl2  }  (4-23) 

Repeated  application  of  the  matrix  inversion  lemma  shows  that 

i r  -  1  -j  ■  (4-24) 
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Again,  under  the  assumption  that  p2°22  »  *>  i-e-.  that  the  power  of  the  second 
jammer  is  well  above  the  noise  and  that  the  projection  of  the  second  jammer’s 
propagation  vector  onto  the  second  basis  vector  is  nonzero 

nah  =  ■  («5| 

This  expression  is,  like  (4-22),  independent  of  the  jammer-to-thermal-noise  ratios 
under  the  assumptions  stated;  it  follows  that  the  ABF  and  CBF  solutions  are 
again  equivalent. 

If  the  methods  and  assumptions  just  described  for  the  single-jammer  and 
dual-jammer  cases  are  extended  to  the  case  of  M  jammers,  then 

Bah  S  </+  E  •  (4-26) 
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That  is,  the  ABF  and  CBF  solutions  become  indistinguishable  for  the  cases  where 
jammer-to-thermal-noise  ratios  in  the  adaptive  processors  remain  large. 


Two  conditions  must  be  satisfied  for  the  jammer-to- thermal-noise  ratio  in 
the  CBF  to  be  large.  In  the  first  place,  of  course,  the  jammer  signal  must  arrive 
at  the  array  elements  with  appreciable  power.  This  will  usually  be  the  case  if  the 
jammer  constitutes  a  threat.  A  second  condition  is  that  the  nullformer  not 
discard  so  much  jammer  energy  that  the  jammer  contribution  is  driven  into  the 
thermal  noise.  In  the  case  of  the  two-element  subarray,  this  difficulty  may  arise 
at  angles  near  the  look  direction,  i.e.,  within  about  10°  of  the  signal  azimuth; 
with  other  nullformers  the  regions  that  seriously  attenuate  jammer  signals  may 
be  different. 

A  third  condition  must  also  be  satisfied  if  the  CBF  weights  are  to  reliably 
approach  optimality.  This  last  condition  relates  to  the  number  of  degrees  of 
freedom  available  versus  the  number  of  jammers  that  must  be  nulled.  In  the 
narrowband  case  one  complex  weight  is  adequate  to  null  a  jammer.  If  there  are 
K  beamformer  channels,  then  the  number  of  jammers  must  be  less  than  or  equal 
to  K- 1  since  one  complex  weight  is  lost  to  the  constraints  in  a  narrowband  Frost 
ABF. 

B.  Further  Simulation  Experiments  with  the  CBF 

Two  additional  simulation  experiments  were  devised  to  assess  CBF 
performance  under  more  difficult  conditions  than  the  single-jammer  scenario  of 
Chapter  3.  Particular  attention  was  given  to  testing  the  optimality  criteria  just 


discussed. 


The  signal/jammer  scenario  for  one  experiment  is  shown  in  Figure  21.  The 
CBF  configuration  is  identical  to  that  used  in  the  experiment  of  Chapter  3,  but 
the  jammer  environment  has  been  altered.  The  jammer  J j  is  a  signal  centered  at 
0.25  ftamp  with  a  bandwidth  of  1%.  The  actual  power  of  Jx  is  2.0,  and  its 
arrival  angle  is  -10*.  The  jammer  */2  is  equal  in  center  frequency  and  bandwidth 
to  Ju  but  it  is  independent  of  Jx.  The  power  of  J2  is  1.5,  and  its  arrival  angle  is 
80*. 

With  the  restricted  number  of  weights  in  the  CBF,  there  were  insufficient 
degrees  of  freedom  to  meet  the  criteria  stated  in  the  previous  section. 
Specifically,  there  were  only  enough  weights  to  null  one  jammer  after  allowing  for 
the  loss  of  degrees  of  freedom  to  the  constraints.  It  was  of  interest  to  examine 
performance  under  these  conditions  of  CBF  saturation. 

Figure  22  shows  the  beampattern  plot  obtained  after  initial  transients  had 
disappeared.  A  relatively  shallow  null  was  placed  on  4  while  there  was  no 
visible  attenuation  of  the  stronger  J  j.  This  behavior  reflected  the  perception  of 
the  jamming  environment  within  the  master  beamformer  where  the  nulling 
circuits  had  attenuated  Jx  and  boosted  J2-  The  actual  solution  reached  by  the 
beamformer  in  this  case  was  a  dynamic  one  in  which  the  power  of  J2  was  reduced 
to  approximate  that  of  Jx.  Additionally,  the  response  to  J2  was  adjusted 
dynamically  to  approximately  cancel  the  contribution  from  Jv  Unfortunately, 
the  solution  was  one  which  mapped  poorly  to  the  slaved  beamformer  because  of 
the  vastly  different  nullformer  influences  on  the  two  jammers. 


ce  Bea 


-  72  - 


A  second  experiment  w as  conducted  in  which  the  same  signal/jammer 
scenario  was  used,  but  the  CBF  was  extended  to  provide  two  additional  degrees 
of  freedom  in  both  the  master  and  the  slave.  Figure  23  shows  the  configuration 
of  this  CBF. 

The  CBF  beampattern  after  convergence  is  illustrated  in  Figure  24. 
Extremely  deep  nulls  have  been  placed  on  both  Jx  and  J2,  and  the  beamformer 
output  consists  almost  exclusively  of  the  desired  signal.  It  is  interesting  to  note 
that  the  null  depth  on  J2  remains  somewhat  greater  than  that  on  Jj.  This 
reflects  the  relatively  greater  jammer- to- thermal-noise*  ratio  of  J2  within  the 
CBF  master  as  a  result  of  nullformer  action. 

C.  Convergence  Time  Constants 

As  was  stated  in  Chapter  2,  convergence  of  the  constrained  LMS  algorithm 
is  governed  by  the  eigenvalues  of  P  Rxx  With  the  aid  of 

P  4  I -C(CTC)-1  CT  (2-16) 

and 

Rxx  —  C  Rrr  CT  +  Rnn  >  (3*2) 

the  matrix  P  Rxx  P  may  be  rewritten  as 

*ThtnpaJ  Boitt  per  tltmtsl  tu  irt  to  —40  iB  relative  to  aait  power. 


p  Rxx  P  =  |/  -  C(CTC)-'  CT\\C  Bkr  Ct  +  Rnn]  P  .  (4-27) 


which  simplifies  to 

P  Rxx  P  =  P  Rnn  P  •  (4-28) 

Once  again  it  is  seen  that  the  role  of  the  desired  signal  in  the  Frost  ABF  is  very 
limited  and  that  the  noise  voltages  govern  convergence.  It  follows  that,  if  an 
ideal  nulling  circuit  were  realizable,  CBF  and  ABF  convergence  rates  would  be 
identical. 

In  practice,  however,  the  nulling  circuits  are  imperfect,  and  the  noise 
voltages  in  the  CBF  master  only  approximate  the  noise  voltages  in  an  equivalent 
ABF.  The  actual  CBF  convergence  rates  are  governed  by  the  noise  voltages  that 
are  present  in  the  master.  These  voltages  may  be  expressed  as 

X  =  ,  (4-29) 

where  A  is  a  matrix  that  describes  the  action  of  the  nulling  circuits  and  Jtf-)  is 
the  vector  of  voltages  from  the  augmented  array  that  is  needed  for  the  CBF. 
Because  the  nulling  circuits  act  to  remove  desired  signal  components  (4-29)  may 
also  be  written  as 

iV*  =  A(S(‘>+  JVW)  =  AJV<*)  .  (4-30) 

The  analog  of  P  Rnn  P  in  the  Frost  ABF  then  becomes 

P  Rnn  P  =  PA  RM)AtP  (4-31) 


for  the  Frost-based  CBF. 
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Nullformer  action  can  either  increase  or  decrease  CBF  convergence  time 
relative  to  a  equivalent  ABF.  An  experiment  was  conducted  to  compare  CBF 
and  ABF  performance  in  the  case  involving  a  single,  sinusoidal  jammer.  The 
ABF  that  was  selected  was  the  two-element  structure  of  Figure  6;  the  CBF  was 
an  equivalent  structure  with  three  elements  and  differencing  circuits.  The  two 
arrays  were  steered  to  broadside  and  constrained  to  provide  an  all-pass  response 
in  that  direction.  Signal  power  was  set  to  zero  to  avoid  the  distracting  influence 
of  signal  cancellation  in  the  ABF.  A  unit-power  jammer  signal  at  a  frequency  of 
0.25  fiamp  was  injected  from  endfire  {0  =  90® ),  and  the  two  systems  were 
allowed  to  adapt  {at  /i  =  1  X  10-2)  until  output  power  was  reduced  to  -20 dB 
relative  to  the  unit  input  power.  The  experiment  was  repeated  to  obtain  the 
number  of  adaptations  to  convergence  at  jammer  angles  ranging  from  endfire  to 
broadside.  The  results  were  plotted  in  Figure  25,  which  shows  the  number  of 
adaptations  to  convergence  versus  jammer  angle  for  both  the  ABF  and  the  CBF. 

Figure  25  shows  that,  depending  upon  jammer  angle,  either  the  ABF  or  the 
CBF  may  provide  the  more  rapid  convergence.  The  CBF  yields  the  better 
performance  in  the  range  from  20*  to  90®  where  the  two-element  array  gain 
exceeds  isotropic.  At  about  20*  the  two  curves  cross,  and  at  angles  nearing 
broadside  the  ABF  provides  superior  performance.  Convergence  of  both  systems 
becomes  very  slow  at  angles  approaching  the  look  direction  where  the 
requirements  of  jammer  nulling  and  the  requirements  of  the  look-direction 
constraints  begin  to  conflict.  The  ABF  will  be  more  effective  in  its  response  to 


Jammer  Angle, 

25.  Comparisons  of  ABF  and  CBF  Convergence  Times  in  the  Single-jammer  Case. 
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jammers  arriving  near  the  signal.  On  the  other  hand,  the  CBF  will  be  somewhat 
more  tolerant  of  small  errors  in  array  steering  because  of  its  sluggish  response  at 
angles  near  the  look  direction.  It  appears  likely  that,  for  small  steering  errors, 
signal  nulling  in  the  CBF  will  be  deferred  indefinitely  by  a  reasonably  dynamic 
electromagnetic  environment  or  by  the  artificial  injection  of  a  small  amount  of 
noise  into  the  system. 

D.  Weight  Quieting  in  the  CBF 

An  important  aspect  of  the  performance  of  an  adaptive  system  is  its 
behavior  near  the  optimal  solution.  Ideally,  the  weights  of  an  adaptive  system 
should  move  to  the  optimal  solution,  remain  at  precisely  that  point  in  weight 
space,  and  displace  only  in  response  to  a  change  in  the  optimal  solution,  i.e.,  in 
response  to  an  actual  change  in  the  environment.  In  practice,  a  compromise 
must  be  reached  between  the  ability  of  an  adaptive  system  to  hold  the  optimal 
solution  in  a  stationary  environment  and  to  track  the  optimal  solution  in  a 
nonstationary  environment.  The  nature  of  the  compromise  is  not  necessarily  the 
same,  however,  from  system  to  system.  This  section  will  compare  the 
performances  of  the  Frost  ABF  and  the  Frost-based  CBF  near  the  optimal 
solution  and  will  show  from  another  point  of  view  that  the  CBF  enjoys  a 
performance  advantage. 

Weight  updating  in  the  Frost  ABF  is  accomplished  in  accordance  with  the 
recursion  relationships  given  in  (2-14)  or  (2-17).  The  same  weight  updating 


technique  is  used  in  both  the  Frost  ABF  and  the  Frost-based  CBF,  but  the  data 
samples  entering  the  computations  differ  significantly  as  the  two  systems 
approach  convergence.  The  contrasts  may  be  seen  by  examining  any  of  the 
weights;  let  the  first  weight  w  t  be  selected. 

Initially  the  Frost  ABF  will  be  considered.  The  difference  between  the 
values  of  w  i  at  two  successive  samples  may  be  written  from  (2-14)  as 

Wj(Hl)-  Wj(Ar)  =  -»y(k)x{(k)  -  £  (wy(*)  -  +  ~- 

(4-32) 

Because  the  constraint  / 1  is,  by  definition,  the  sum  of  the  weights  in  the  first 
column  of  the  beamformer,  (4-32)  may  be  simplified  to 

Wi(*+1)-  w,(*)  =  -ny{k)  [*,(*)  -  jr  £  xy(*))  . 

(4-33) 

Let  it  be  assumed  that  at  the  kth  sample  instant  the  adaptive  system  has 
reached  a  set  of  weights  that  closely  approximates  the  optimal  solution.  Then 

w,(*)  =  w,#  (4-34) 

Additionally,  if  the  constraints  are  set  to  provide  unity  gain  and  an  all-pass 
response  in  the  look  direction,  the  system  output  is 

y(k)  =  a{k)  +  nr{k)  ,  (4-35) 

where  s(k)  is  the  signal  at  the  kth  sample  time  at  any  of  the  taps  in  the  first 
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column  of  the  beamformer  and  nr(k)  is  the  residual  noise  that  exists  at  the 
beamformer  output  with  optimal  weights. 

The  tap  voltages  in  the  first  beamformer  column  may  also  be  stated  in  terms 
of  signal  and  noise  components: 

x,(fc)  =  «(*)  +  n,(fc)  ,  i  =1,  2,  K  (4-36) 

Note  that  the  cophased  signal  is  the  same  at  taps  1  through  K.  This  allows  (4- 
33)  to  be  simplified  to 

[wi(*+  1)  -  w^*)]^  =  -/i(*(*)  -I-  nr(k)\  (n,  -  ~  £  nyj  (4-37) 

A  /-l 

The  CBF  weight  recursion  will  now  be  considered  under  like  circumstances. 
Again  assume  operation  near  W*.  Further  assume  that  the  nulling-circuit 
performance  is  nearly  ideal  so  that  differences  in  the  noise  voltages  between  the 
CBF  and  a  comparable  ABF  may  be  neglected.  Then  the  CBF  tap  voltages  are 

*,(*)  =  «,•(*)  ,  «  =  1, 2 . K  .  (4-38) 

The  weight  difference  is  obtained  from  (4-33)  as 

lwi(*+l)-w  i{k)]CBF  =  -/*[»r(*))  ("i  “  -h  E  n/l  (4-39) 

K  1 

The  contrast  between  ABF  and  CBF  weight  behavior  near  convergence  may 
be  seen  by  comparing  (4-37)  and  (4-39);  the  difference  lies  in  the  inclusion  of  the 
signal  term  in  the  ABF  update  expression. 


The  significance  of  this  contrast  is  illustrated  in  Figure  26,  which  shows  a 
cut  along  the  w  |  axis  through  the  error  surfaces  of  the  two  systems.  Under  the 
assumption  of  negligibly  different  noise  voltages  in  the  ABF  and  CBF,  the  shapes 
of  the  two  error  surfaces  are  identical.  The  ABF  minimum,  however,  is 
substantially  elevated  above  the  CBF  minimum.  If  a  useful  signal  is  to  be 
delivered  by  the  adaptive  system,  output  signal-to-noise  ratios  in  excess  of  10  dB 
are  ordinarily  required,  and  ratios  in  excess  of  20  dB  are  certainly  desirable. 
These  ratios  imply  that  the  signal  levels  will  typically  be  three  to  ten  times  the 
noise  levels  as  the  ABF  nears  convergence.  The  typical  weight  adjustment  in  the 
ABF  will  thus  be  dominated  by  signal  effects  as  the  adaptive  system  nears 
convergence.  This  becomes  especially  important  when  it  is  realized  that  the 
signal-times-noise  terms  amount  to  short-term  estimates  of  the  correlation 
between  the  signal  and  the  noise.  In  other  words,  the  signal  influence  in  the  ABF 
should  not  be  regarded  as  simply  a  random  weight  adjustment,  but  as  a 
purposeful  step  toward  effecting  signal  cancellation.  Because  the  CBF  weight 
adjustments  are  not  influenced  by  the  signal,  the  CBF  weights  exhibit  less 
variation  about  the  optimal  solution  than  the  weights  of  the  comparable  ABF. 
Further,  the  CBF  weight  adjustment  process  cannot  cause  signal  cancellation. 


V.  SUMMARY  AND  CONCLUSIONS 


It  has  been  shown  that  interaction  between  the  desired  signal  and 
interference  signals  can  lead  to  partial  or  total  cancellation  of  the  desired  signal 
within  an  adaptive  beamformer  (AGF).  The  signal  cancellation  phenomenon  may 
arise  spontaneously  in  a  multipath  environment,  or  it  may  be  deliberately  or 
accidentally  induced  by  intentional  or  unintentional  jamming.  High  adaptation 
rates  increase  the  susceptibility  to  signal  cancellation  in  beamformers  that  have 
not  been  designed  to  avoid  the  effect. 

A  composite  beamformer  (CBF)  has  been  described  that  avoids  signal 
cancellation.  The  CBF  makes  use  of  both  an  adaptive  beamformer  and  a 
nonadaptive,  slaved  beamformer.  The  inputs  to  the  adaptive  beamformer  are 
preprocessed  to  exclude  the  desired  signal,  but  retain  representation  of  the 
interference  signals.  By  restricting  the  adaptive  process  in  this  fashion,  a  solution 
(i.e.,  a  set  of  beamformer  weights)  is  generated  that  influences  only  the 
interference  signals.  Weights  from  the  adaptive  processor  are  copied  to  the 
slaved  processor,  which  is  connected  directly  to  the  antenna  elements  and  which 
produces  the  useful  output  of  the  array.  The  CBF  delivers  much  better 
performance  with  regard  to  signal  cancellation  than  a  more  conventional  ABF. 
Simulation  data  have  been  presented  to  show  the  improvement  that  can  be 
obtained. 

The  performance  of  a  Frost-based  CBF  has  been  compared  with  that  of  a 
Frost  ABF.  If  a  narrowband  system  is  assumed,  it  has  been  shown  that  the 


9.  —  * 


.  'JL TmA^m 


»  *  w  * 


-  84  - 


Frost-based  CBF  will  tend  toward  the  optimal  (i.eM  the  Wiener)  solution 
provided  1)  the  preprocessor  admits  energy  from  each  jammer,  2)  energy  from 
each  jammer  dominates  the  thermal  noise,  and  3)  there  are  sufficient  degrees  of 
freedom  to  null  each  jammer.  The  properly  designed  CBF  is  thus  equivalent  to 
the  comparable  ABF  in  terms  of  the  ideal  solution.  The  CBF  will  not,  however, 
exhibit  signal- induced  weight  variations.  Performance  will  therefore  be  improved 
in  those  cases  where  signal  cancellation  occurs  in  the  ABF.  Convergence  toward 
the  optimal  solution  will,  in  general,  be  different  for  a  CBF  and  an  equivalent 
ABF  that  are  operating  in  the  same  environment.  The  convergence-rate 
differences  are  governed  by  the  influence  of  the  preprocessor.  In  some  cases  the 
CBF  will  attenuate  a  given  interferor  more  rapidly  than  the  ABF,  while  in  other 
cases  the  reverse  is  true.  If  the  preprocessor  is  designed  to  match  the  basic  array 
element  pattern  except  for  a  sharp  null  in  the  look  direction,  the  convergence 
rates  for  the  CBF  and  the  ABF  are  essentially  the  same.  An  exception  occurs  for 
a  jammer  incident  from  very  near  the  look  direction.  The  CBF  responds  even 
more  sluggishly  than  the  ABF  to  such  an  interferor.  It  is  also  true,  however,  that 
the  CBF  is  less  sensitive  to  small  steering  errors  than  the  ABF. 

The  studies  of  signal  cancellation  and  the  CBF  have  thus  far  suggested 
several  areas  for  further  investigation.  The  signal  cancellation  phenomenon  itself 
is  not  so  well  characterized  as  might  be  desired.  An  improved  characterization  of 
the  phenomenon  may  provide  a  better  means  for  analyzing  the  importance  of 
signal  cancellation  in  various  systems  and  environments.  Additionally,  this  work 
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may  reveal  new  ways  to  avoid  the  effect  in  adaptive  arrays. 

Perhaps  the  most  promising  area  for  further  CBF  improvement  is  in  the 
preprocessor  that  nulls  the  target  signal.  Work  thus  far  has  centered  on  the 
element-differencing  scheme  because  of  its  simplicity.  Other  preprocessor 
configurations,  though  more  complex  from  a  hardware  viewpoint,  are  capable  of 
yielding  subarray  patterns  that  provide  a  better  match  between  the  patterns  seen 
by  the  adaptive  beamformer  and  the  slaved  beamformer.  Improved  matching 
will  assure  that  the  CBF  performance  (both  dynamic  and  asymptotic)  is 
approximately  that  of  the  equivalent  ABF,  except  for  the  superior  CBF 
performance  with  regard  to  signal  cancellation.  Steyskal  [21]  has  recently 
described  a  procedure  for  optimal  pattern  matching  subject  to  null  constraints. 
This  procedure  appears  to  be  well  suited  for  preprocessor  design. 

The  general  procedure  of  adapting  in  the  absence  of  signal  can  sometimes  be 
applied  without  resorting  to  the  CBF.  Systems  that  make  use  of  low-duty-cycle 
or  frequency-hopped  signals  can  employ  other  signal-exclusion  schemes  that  arise 
naturally  out  of  the  signaling  format.  Work  is  in  progress  at  Stanford  on 
cancellation-free  adaptation  schemes  for  systems  in  these  classes.  These  schemes 
will  trade  timing  constraints  against  the  spatial-nulling  method  of  the  CBF  and 
will  almost  certainly  result  in  less  complex  adaptation  schemes  for  systems  that 
do  not  continuously  occupy  a  channel. 
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APPENDIX  A 

ANALYSIS  OF  THE  DYNAMIC  SOLUTION 
FOR  THE  SINUSOIDAL  CASE 

Chapter  2  discussed  the  dynamic  solution  that  is  tracked  by  an  ABF  when 
large  values  of  the  adaptation  parameter  ;s  are  used.  A  computer  simulation 
formed  the  basis  of  that  discussion,  and  certain  claims  were  made  without  a 
supporting  analysis.  In  this  appendix  the  dynamic  error-surface  minimum  is 
examined  analytically,  and  a  mathematical  model  is  developed  that  supports  the 
discussion  in  Chapter  2. 

In  Chapter  2  the  Frost  ABF  under  study  was  mapped  into  an  equivalent 
adaptive  noise  canceller  (ANC)  in  order  to  reduce  the  problem  dimensionality  and 
simplify  the  analysis.  This  same  mapping  will  be  used  here;  Figure  A-l  reviews 
the  ANC  structure  that  will  be  analyzed. 

The  ANC  structure  assumes  separate  availability  of  the  desired  signal  and 
the  single  jammer  signal,  both  of  which  are  sinusoids  in  the  case  under 
consideration.  The  desired  signal  is  taken  directly  to  the  output  summing 
junction  in  conformance  with  a  unit-gain,  all-pass  constraint  in  the  look 
direction.  The  jammer  signal  is  processed  through  a  two-channel  (in-phase  and 
quadrature)  filter  with  one  weight  per  channel  and  is  subtracted  from  the  desired 
signal  at  the  output  summing  junction.  The  .ANC  output  is  also  the  error  signal 
that  is  fed  back  to  control  adaptation  of  the  filter  weights.  In  the  notation 
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Figure  A-l.  Equivalent  Adaptive-Noise-Canceller  Representation. 
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indicated  in  the  figure,  the  error  is 

t{x)  =  d(x)  -  K*i(*)  +  w2*2(z)l  .  (A-l) 

The  task  of  the  adaptation  algorithm  is  to  minimize  some  function  of  e(z) 
by  selecting  appropriate  values  for  Wj  and  w2.  Under  the  assumptions  of 
stationarity  and  slow  adaptation,  the  effectiveness  of  an  adaptation  algorithm  is 
often  established  by  comparing  performance  with  a  minimum-mean-square-error 
solution  that  is  derived  analytically.  It  has  been  established  that  the  Widrow- 
Hoff  LMS  algorithm  delivers  asymptotically  optimal  performance  in  the 
stationary  case,  i.e.,  the  LMS  solution  tends  to  the  Wiener  solution  as  p  tends  to 
zero. 

Here,  however,  interest  centers  upon  an  ABF  that  is  adapting  rapidly  and  is 
therefore  operating  on  the  basis  of  short-term  estimate*  of  the  environment.  At 
high  adaptation  rates  it  is  unrealistic  to  expect  close  conformance  with  the 
Wiener  solution,  which  is  derived  under  a  very  different  view  of  the  data. 
Instead,  an  objective  function  is  needed  that  more  accurately  reflects  the  use  of  a 
finite  time  window  and  allows  the  solution  to  evolve  over  time.  An  alternative 
objective  function  (  will  be  taken  as  the  average  square  error  over  the  interval  A: 

'at  a)  k  —  jr^V)  *  •  (a-2) 

This  objective  function  has  the  desirable  property  of  extending  to  mean-square 
error  as  t  and  A  increase  without  bound.  Nebat,  et  al,  [1]  have  made  use  of  a 
similar  finite-time  integration  in  exploring  the  effect  of  record  length  on  the 
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correlation  of  complex  exponentials. 
Combining  (A-l)  and  (A-2)  gives 


t(*A)  =  X  f*J2(x)dx  -  ft[Ad(z)zl(z)dz  -  ^  ft'Ad(z)z2(z)dx 


2 to.  .t 


2w2  ,* 


wf  . t  2w.w o  f  w« 

+  “  +  ~ ~  J[-4  *1<*W*)*  +  -3-  J,_4  *!(*)* 


(A-3) 


To  provide  a  more  compact  notation,  let 


1  r* 


rij(tA)  4  X  /|-A  I,  y  =  1,  2 


Pit*  A)  4  J[*A  d(*K(x)d»  »  =  1,  2 


(A-4) 


Equation  (A-3)  may  now  be  rewritten  as 


£(/,A)  =  /{  A  d2(z)dz  -  2 to,  p,(t,A)  -  2to2  p2(*»A) 

+  wf  r„(f,A)  +  2u>,w2  rl2(f,A)  +  w|  r^/.A)  £A-5) 


* 

Minimization  of  £(f,A)  is  accomplished  by  taking  the  gradient  with  respect 
to  the  weights,  setting  the  gradient  to  zero,  and  solving  for  w ,  and  w  2.  The  first 
two  steps  of  this  procedure  yield 

[-2  Pi(tA)  +  2 to,  r„(/,A)  +  2to2  rl2(t, A)  1  |o  ] 

,A)  =  l2  n(tA)  +  2Wi  rK{tA)  +  2Wz  f22{tA)  I  -  [0  I  (A-6) 


The  solution  of  (A-6)  may  be  written  as 


rvro 


- ’ o . »  «V 


cm 
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>n(*.A)  ^.A) 

-I 

Pi(*,A) 

»2l(<.A)  r^A) 

,P2(<»A) 

(A-7) 


where  wlm  and  w2m  are  the  values  of  the  weights  that  minimize  the  objective 
function  during  a  given  time  window  and  where  existence  of  the  matrix  inverse  is 
assumed.  The  equivalent  matrix  equation  for  (A-7)  is 


Wm  =  R'HtA)  P{‘A) 


(A-8) 


Equations  (A-7)  and  (A-8)  are  finite-time-window  versions  of  the  familiar 
Wiener-Hopf  equation. 

Since  the  dimension  of  (A-7)  is  only  two,  it  is  simple  to  expand  the  matrix 
inverse  and  obtain  expressions  for  the  weights: 

jctmM )  {,a(,’A)  '  r,a“'A)  r*tA) } 

]  |-r,,(<>A)p1(l,A)+  r„(<,A)  P2(*,A)j  .  (A-0) 

Equations  (A-0)  describe  the  optimal  weights  for  the  two-weight  adaptive  noise 
canceller  at  time  t  when  averaging  is  over  a  window  of  duration  A. 

In  order  to  obtain  a  direct  comparison  with  the  sinusoidal  case  considered  in 
Chapter  2,  let 


d(x)  as  COSU|X 
xt(x)  —  cos  (w2*+  <t>) 
z2(x)  =  sin  (w2*+  <j>) 


(A- 10) 
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These  signal  definitions  lead  to  the  following  values  for  components  of  the 
j?(f  ,A)  matrix  and  the  P(t,A)  vector: 


r^tA) 

Pi(*»A) 


-  i+i  1  e“[*th + 

=  r2i(‘A)  =  j  s*nc|“““  J  **•  *  +  ] 

-  »  - *  ,to,!hr)cos[!!h,  +  *+ 

1  .  f(“i-W2)A  1  f,  ,  .  (w,-u2)A 

=  7  —f-5—  I  “T 


7  sinc[ 


pd*A) 


=  -  —  sine 
2 


(Wi-ujaJA 
2  x 

(wj+o^A 
2x 

((tJj-id 2)A 

2* 


I(u7|-w2)/  -  <f>  + 


cos|(W|+fa>2)f  +  <l>  + 


(wj+  w2)A 


sin  |(u/i-w2)f  -  <f>  + 


(w,-w2)A 


I  (w.+  w2)A 

sine  . .  .  sm  (wj+  u2)t  +  <t> 

&  I/K 


(w,+  u/2)A 


(A-ll) 


Equations  (A-9)  and  (A-ll)  provide  a  complete  specification  of  ideal  weight 
behavior  in  the  ANC  problem  being  considered.  It  would  be  possible  at  this 
point  to  utilize  these  equations  to  write  expressions  for  wlm  and  w2m.  It  is 
apparent,  however,  that  the  combined  complexities  of  (A-9)  and  (A-ll)  would 
lead  to  a  rather  unwieldly  pair  of  equations  for  the  weights.  A  more  satisfactory 
approach  is  to  continue  with  the  existing  equations  and  examine  some  particular 
aspects  of  weight  behavior  that  bear  upon  the  discussion  in  Chapter  2. 

Initial  attention  will  be  directed  toward  Equations  (A>11).  Sine  functions 
appear  prominently  in  the  expressions  for  all  of  the  components  of  the  R  (f  ,A) 
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matrix  and  the  P(t, A)  vector.  These  sine  functions  govern  the  peak  amplitudes 
of  most  terms  in  the  correlation  expressions  and  thus  play  a  key  role  in 
determining  weight  behavior.  The  arguments  of  the  sine  functions  involve  only 
the  time  window  A  and  some  function  of  u>l  and  u2.  This  latter  quantity  is  fixed 
for  a  given  problem  statement,  but  the  equivalent  time  window  can  be  governed 
in  an  adaptive  system  by  choice  of  the  adaptation  parameter  #i.  The  observation 
that  weight  amplitudes  can  be  influenced  by  the  choice  of  adaptation  parameter 
harmonizes  well  with  the  behavior  observed  in  the  experiments  of  Chapter  2. 

A  better  understanding  of  the  relative  importance  of  the  terms  in  (A-ll)  can 
be  obtained  from  a  closer  study  of  the  sine  functions.  Figure  A-2  is  a  plot  of  sine 
function  magnitude  as  a  function  of  the  argument.  An  upper  bound  for  the  sine 
function  is  also  plotted  to  reflect  the  inverse  dependency  on  argument  that 
dominates  the  sine  function  behavior.  It  is  clear  from  the  figure  that  a  small 
argument  is  necessary  if  the  sine  function  is  to  contribute  significant  amplitude  to 
a  term.  Conversely,  it  is  easy  to  see  that,  unless  Wj  and  w2  are  equal,  large  values 
of  A  will  lead  to  large  arguments  and  vanishing  amplitude  contributions  from  the 
terms  involving  sine  functions,  i.e., 

lim  r„(f,A)  =  \ 

lim  r12(f,A)  =  lim  r21(f,A)  =  0 

A -♦oo  A-*oo 

lim  r^f.A)  =  \ 

A-*oo  4 
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iim  Pi(*,A)  =  0 

A-*oo 

Iim  p2 (t,A)  =  0  .  (A- 12) 

A-»oo 

These  are  the  same  values  of  the  correlation  functions  that  would  have  been 
obtained  from  the  traditional  calculations  involving  expected  values.  The  values 
in  (A-12)  may  be  used  in  (A-9)  to  compute  the  Wiener  weights: 

w  =  0 

w  j*  =  0  .  (A- 13) 

The  vector  components  pt(f , A)  and  p2(f  ,A)  each  involve  sine  functions  with 
arguments  containing  the  sum  u>x  +  u2  and  difference  u>x  -  w2  of  the  signal  and 
jammer  frequencies.  Since  ux  and  u2  will  be  taken  as  non-negative  numbers  in 
this  context, 

wi  +  w2  ^  ui  ~  u2  (A-14) 

for  all  allowed  choices  of  and  u2.  When  appreciable  levels  of  cancellation  are 
occurring,  the  argument  A(w1-w2)/2j  will  be  small  so  that  the  leading  terms  in 
Pi(f,A)  and  p^tA)  will  have  magnitudes  approaching  unity.  The  argument 
A((jj+  w2)/2x  can  be  no  smaller  than  A(w1-u;2)/2t,  and  can  be  much  larger.  In 
the  case  considered  in  Chapter  2,  for  example,  wt  +  ui2  was  approximately  SO 
times  larger  than  ux  -  w2.  This  difference  was  sufficiently  large  to  assure  that 
difference-frequency  terms  dominated  weight  behavior  over  the  entire 
adaptation-parameter  range  that  yielded  measurable  deviation  from  the  Wiener 
solution.  Sum-frequency  influence  was  visible  in  the  weight  loci  for  some  values 
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of  the  adaptation  parameter,  but  never  caused  deviations  of  more  than  a  few 
percent  from  the  basically  circular  loci. 

It  is  useful  to  look  more  closely  at  weight  behavior  when  difference-frequency 
terms  are  dominant.  Assume  that  and  w2  are  nearly  equal  and  are  sufficiently 
different  from  zero  so  that 


(w,-w2)A  _  „ 

s—  < 0  5 


(A-15) 


and 
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(A*  18) 


Now,  from  (A-9) 


w =  sine 


w2m  =  atnc 


(wj— u/2)  A 


(Wl-U>2)A 


j  cosl( 


(W|— w2)f  -  <f>  + 


Jsin|( 


(w,-w2)t  -  ^  + 


(w,-w2)A 


(co,-w2)A 


(A-19) 


Equations  (A-19)  describe  a  phasor  rotating  in  the  w1  -  w2  plane  at  the 
difference  frequency  (wj-w2).  The  magnitude  of  the  phasor  is  proportional  to  the 
magnitude  of  the  cross-correlation  functions  pt(f  ,A)  and  p2(f  ,A). 

For  a  given  time  window  A  and  a  specified  time  t ,  the  phasor  given  by 
Equations  (A-19)  defines  the  weight  pair  that  minimizes  the^ANC  output. 
Equivalently,  the  phasor  marks  the  projection  of  the  dynamic  error  surface 
minimum  onto  the  Wj  -  w2  plane.  As  time  evolves,  the  tip  of  the  phasor  traces 
out  the  circular  locus  of  the  projection  of  the  minimum.  This  is  the  behavior 
that  was  exhibited  in  Figure  12  of  Chapter  2. 
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APPENDIX  B 

AN  ALTERNATIVE  EXPRESSION  FOR  THE 
OPTIMAL  WEIGHTS  IN  A  FROST  ABF 


Frost  (1,2]  derived  an  expression  for  the  optimal  weight  vector  in  a  hard* 


constrained  adaptive  beamformer: 


WUr  =  C\CTRi)(C\-'  l  ,  (B-l) 

where  R %%  =  £(.XXr]  is  the  correlation  matrix  for  the  tap  voltages, 
L  =  l/ii  /2t  — *  f  j\T  **  a  constraint  vector,  and 


b  a  constraint  matrix.  In  thb  expression  for  the  dependency  on  the  signal 

environment  b  given  in  terms  of  the  tap  voltages,  i.e.,  in  term  of  Rxx>  which 
may  be  estimated  from  X. 

In  thb  appendix  it  will  be  shown  that  an  equivalent  formulation  for  the 
optimal  weight  vector  may  be  given  in  terms  of  the  vector  N  where 

N  =  X-S  (B-3) 

b  the  vector  of  tap  nobe  voltages  that  b  obtained  by  deleting  the  desired-signal 
component  5.  In  particular,  it  will  be  shown  that  the  equivalent  expression  b 

H 'itr  =  l  ,  (B-4) 

where  RNN  ~  K[:V;V7"|  is  the  correlation  matrix  for  the  tap  noise  voltages. 

The  balance  of  thb  appendix  b  broken  into  two  sections.  In  the  first  section 
the  validity  of  (B-4)  will  be  demonstrated  for  the  narrowband  case.  The  second 
section  will  then  examine  the  somewhat  more  complicated  wideband  case. 

A.  The  Narrowband  Case 

In  the  narrowband  case  it  b  assumed  that  signal  bandwidths  are  small 
enough  to  allow  incident  signab  to  be  characterized  as  sinusoids  with  amplitude 
and  phase  variations  that  are  slow  with  respect  to  the  carrier  frequency.  Under 
these  circumstances  a  single  complex  weight  b  adequate  to  adjust  the  gain  and 
phase  of  each  beamformer  channel.  That  b,  the  tapped-delay-line  filter 
constituting  each  beamformer  channel  b  replaced  by  an  in-phase  weight  and  a 


quadrature  weight.  For  the  narrowband  Frost  ABF,  the  constraint  vector 
becomes  simply 

l  =  ft  ,  (B-5) 

with  / 1  a  complex  number  giving  the  in-phase  constraint  as  its  real  part  and  the 
quadrature  constraint  as  its  imaginary  part.  Additionally,  the  constraint  matrix 
takes  a  simple  form: 

C  =  [1  1  •  •  •  l]r  4  1  (B-6) 

The  correlation  matrix  Rxx  must  be  decomposed  into  its  desired-signal 
component  and  its  noise  component  in  order  to  proceed  further.  This  may  be 
accomplished  by  rewriting  (B-3)  in  the  form 

X  =  N  +  S  (B-7) 

and  then  making  use  of  the  cophased  nature  of  the  signal  to  obtain 

X  =  N  +  v7  i  ,  (B-8) 

where  the  signal  amplitude  is  taken  as  \T»  .  Now,  (B-8)  may  be  used  to  obtain 

Rxx  =  E[XXT ]  =  £[AWr]+  aiir 

=  Rm  +  *1  lr  >  (B-9) 

where  Frost’s  assumption  [1]  of  a  lack  of  correlation  between  signal  and  noise 
voltages  has  been  used  to  simplify  the  expression  for  Rxx- 

It  will  be  useful  at  this  point  to  focus  attention  on  the  factor  Rjix  C  that 


appears  twice  in  (B-l)  and  to  make  use  of  an  observation  by  Baird  and  Zahm  [3] 
regarding  the  utility  of  the  matrix  inversion  lemma  in  a  closely  related  problem. 
The  lemma  may  be  drawn  from  (4]  in  the  following  form: 


If  a  nonsingular  matrix  H  of  order  n  can  be  written  in  the  form 


H  =  A  +  BDE  ,  (B-10) 

where  B ,  D,  and  E  are  n  X  p,  pX  p,  and  p  X  n,  respectively,  and  A  and 
D  are  nonsingular,  then 

H  l  =  A~l  -  A~lB(D~i+  EA~XB)~X  EA~l  .  (B-ll) 

The  lemma  given  in  (B-ll)  may  be  used  to  invert  Rxx  in 

RAC  =  I  Rm  +  .nr]-‘  1  ■  (B-12) 

As  a  preliminary  step  it  is  noted  that  Rfjfj,  which  is  analogous  to  the  A  matrix 
in  (B-10),  is  invertible  because  it  is  assumed  to  be  positive  definite.  The  vectors 
«1  and  lr,  which  correspond  to  B  and  E  in  (10),  are  K  x  1  and  lx  K, 
respectively.  The  D  matrix  in  (B-10)  degenerates  to  the  scalar  1  in  (B-12)  and  is 
clearly  invertible.  With  the  applicability  of  the  inversion  lemma  established,  (B- 
12)  may  be  written  as 


Rx'xC 


9  Rah  1  lr  Rm 
(I  +  *  lr  Rnn  I) 


1 


(B-13) 


Expanding  gives 
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B.  The  Wideband  Case 

The  beamformer  in  the  wideband  case  is  assumed  to  contain  K  filters,  each 
with  J  taps.  The  AV-element  vector  X  of  tap  voltages  may  again  be  written  as 

X  =  N  +  S  ,  (B-18) 

where  N  is  the  vector  of  tap  voltages  attributable  to  noise  and  interference  and 
S  is  the  desired-signal  vector.  In  general,  the  components  of  N  will  each  be 
different  due  to  independent  noise  in  each  channel  and  due  to  the  fact  that 
interference  signals  will  typically  not  be  aligned  across  the  array.  The  signal,  on 


the  other  hand,  is  aligned  across  the  array  and  is  therefore  identical  down  any 
column  of  taps  within  the  beamformer.  The  special  structure  of  S  may  be 
expressed  by  the  relationship 


S  =  CR  ,  (R-19) 

where  C  is  defined  by  (B-2)  and 

R  4  [«(*)  s(k-r)  »(k-2r)  •  •  •  s(l k-  (J-l)r)  ]r  (B-20) 

is  the  vector  of  desired-signal  voltages  at  the  taps  of  any  beamformer  delay  line. 
Combining  (B-18)  and  (B-19)  yields 

X  =  N+  CR  .  (B-21) 

The  correlation  matrix  for  the  wide  band  case  may  be  computed  based  upon 
(B-21): 

Rxx  =  £[(W  +  Ci?)(^r+  RtCt)  )  (B-22) 

Rxx  =  E\NNt\  +  CE\RRt)Ct  (B-23) 

The  transition  from  (B-22)  to  (B-23)  has  made  use  of  the  assumed  lack  of 
correlation  between  signal  and  noise.  Alternative  forms  for  (B-23)  are 

Rxx  —  Rnn  +  c  RrrCT  (B-24) 

and 

Rxx  1=5  RffN  +  Rss  •  (B-25) 

It  should  be  noted  that  RN N  and  Rss  are  dimensioned  KJ  x  KJ,  while  Rrr  is 


dimensioned  J  x  J.  The  difference  in  dimension  between  Rss  and  Rrr  reflects, 
of  course,  the  redundancy  in  Rss- 


The  matrix  inversion  lemma  that  was  used  in  the  narrowband  case  will 
again  be  used  to  invert  the  matrix  Rxx  as  expressed  in  (B-24): 

Rjoc  —  Rnn  ~  Rss^[Rrr  +  Cr  ^jvArC'J'1  CtRn'n  (B-26) 

Both  Rf/ff  and  Rgg  are  assumed  to  be  positive  definite  and  therefore  invertible. 

The  relationship  given  in  (B-26)  will  now  be  used  to  structure  components  of 
the  expression  for  W^gF  given  in  (B-l).  The  leading  factor  is  easily  constructed 
by  post  multiplying  (B-26)  by  C: 

R&C  =  R^C  -  Rjjhcml  +  CrKjjJ,C|-'  CtRn'„C  .  (B-27) 

The  central  factor  is  obtained  in  two  steps.  First,  (B-27)  is  premultiplied  by 
CT: 

CtR&C  =  ICrfl*J,C]  +  |-Crfl^C)  . 

(B-28) 

The  matrix  inversion  lemma  is  then  applied  once  again  to  obtain: 

I CtRx!cC]-'  =  |Crfl4C|-'  +  |<7^C|  {  (fi*i  +  CtR^C\ 

-  |CrB^C|[Cr/?Jv]lrC|-'  [<7rfl^C|}-> 
[C'R^CWCTR^C]-' 

\CtR£C\-'  =  +  Rss  .  (B-20) 


Equations  (B-27)  and  (B-29)  may  now  be  used  to  synthesize 
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Rj^cC\CT Rj^[C\~l  =  {RfihC  -  Rf^ClR^h  +  ^tRnnC]~1  CtR^C} 

(I  CTRti\iC\  +  Rrr)  (B-30) 

Performing  the  indicated  multiplication  on  the  righthand  side  of  (B-30)  and 
simplifying  yields 

RxxC\CT Rj&C\~l  =  Rnn  C[ CTRf^C]~l  -4-  RjilfC  Rrr 

-  RnnC[Rb)i  +  CTRi/lfCYxh 

-  +  CtRn\iCYx  CTRmc  Rrr  (B-31) 

After  it  is  recognized  that  the  J  x  J  identity  matrix  //  may  be  factored  as 

h  =  Rih  Rrb  i  (B-32) 


(B-31)  may  be  rewritten  as 

RxxC\CT Rj&C\~x  =  Rn\iC\Ct Rik^Y1  +  RmC  Rrr 
-  RnMRA  +  Crf?^C}->  [R&  +  CtRJnC\  Rrr  .  (B-33) 

This  last  equation  readily  simplifies  to 

RxxC\CT Rj&C\~x  =  R*kC[CTRs'NCY'  (B-34) 


Postmultiplying  (B-34)  by  the  constraint  vector  /  leads  to  the  desired  result: 

WXbf  =  RxxC\CT Rxx~lCYx  X  =  RnnC\CTRnnCYx  X 
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APPENDIX  C 


AN  ANALYSIS  OF  SIGNAL  INFLUENCE  BASED 
UPON  THE  GENERALIZED  SIDELOBE  CANCELLER 

Another  vehicle  for  analyzing  the  role  of  the  desired  or  look-direction  signal 
b  provided  by  the  Griffiths-Jim  generalized  sidelobe  canceller  [1,2],  a  form  of 
which  b  diagrammed  in  Figure  C-l.  In  thb  beamformer  there  b  a  fixed  processor 
(shown  at  the  right  side  of  the  figure)  that  operates  directly  upon  the  signab 
from  the  array  and  implements  the  look-direction  constraints  as  stated  by 
{/'.},  i  si,  2,  ...,  N.  Note  that  the  look  direction  b  perpendicular  to  the  array 
in  the  case  shown;  steering  delays  are  needed  for  other  look-direction  choices. 
There  b  also  a  second,  adaptive  processor  that  operates  upon  the  difference 
signab  between  adjacent  elements  and  therefore  b  not  exposed  to  look  direction 
signab  through  the  signal  vector  X! .*  The  output  of  the  adaptive  processor  b 
subtracted  from  the  output  of  the  fixed  processor  to  obtain  a  useful  signal.  Thb 
signal  also  constitutes  the  error  signal  that  b  used  in  an  LMS  update  algorithm 
as  opposed  the  more  complicated  constrained  LMS  algorithm  that  would 
otherwise  be  needed  to  implement  hard  constraints.  It  has  been  shown  [2],  that 
the  Frost  ABF  and  the  generalized  sidelobe  canceller  generate  equivalent 
solutions  when  the  sidelobe  canceller  takes  the  particular  form  shown  in 

*Tkere  b,  kowrrtr,  ex  poet  re  to  tit  look-directioa  eipal  tkroifi  tit  error  tiptl  that  b  fed  bock  to  direct  adap¬ 
tation,  tie  (eaerafutd  tidtiobe  caaeeUer  b  •abject  to  ripal  caieellatioa. 


«*?**%* 


Figure 


Figure  C-l. 


In  terms  of  the  difference  signal  vector  X *  and  the  signal  d,  the  Wiener 
weight  vector  for  the  generalized  sidelobe  canceller  may  be  written  as 

W£sc  =  iX*  (C-l) 

where 

x'  41*;,*; . 4/ 

4  E\X‘X'T) 

Pjc-i  4  B\dX‘\  .  (C-2) 

A  representative  element  x- ,  l  <  i  <  A-l,  of  the  difference  signal  vector 
X'  is  obtained  from  the  received  signals  according  to 

xi  =  1  •  (C-3) 

(The  remaining  elements  of  X '  were  derived  in  the  same  manner  at  an  earlier 
point  in  time.)  Each  received  signal  is  taken  to  consist  of  a  signal  component  plus 
noise  components.  Because  the  array  is  assumed  to  consist  of  identical  elements 
and  to  be  perfectly  steered,  the  signal  components  are  identical  across  the  array. 
Therefore 

xi  =  =  («+  n.  )  -  («+  n,.j)  =  n,  -  nM  .  (C-4) 

It  is  evident  from  the  lack  of  signal  components  in  the  difference  signal  vector 
that  the  factor  R^x*  can  bring  no  signal  influence  to  the  Wiener  solution. 


The  output  d  of  the  fixed  beamformer  in  the  generalized  sidelobe  canceller  is 
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a 


m 


given  by 


K  i 

E 

.-i  K 


—  "b  E  (*+  **»)  =  *  +  \  E  n%  (C“5) 

A  »-i  K  •-! 


The  crosscorrelation  may  be  calculated  from  X '  and  d  as  follows: 

P'iX  4  E[dX>) 

=  *[(.  +  -  E  «.)*') 

K  i— 1 

=  £(«X'1  +  £?((-:  E 

A  .-i 

=  E[(j-  E  «,)*'!  • 

*  1-1 


(C-6) 


It  is  now  seen  that  the  second  factor  of  Equation  (C-l)  b  also  free  of  signal 
influence  and  therefore  that  the  Wiener  solution  for  a  hard-constrained  ABF  b 
uninfluenced  by  the  desired  signal. 
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APPENDIX  D 


A  COMPOSITE  BEAMFORMER  BASED  UPON 
WIDROW’S  PILOT-SIGNAL  ALGORITHM 

The  composite-beamformer  discussions  in  the  body  of  this  report  centered 
upon  the  use  of  Frost's  algorithm  for  weight  adaptation.  Frost’s  algorithm  is 
designed  to  operate  in  a  signal-aligned  environment,  i.e.,  the  array  must  be 
steered  either  mechanically  or  electrically  so  that  signals  from  the  look  direction 
are  cophased  across  the  array.  Once  the  array  has  been  properly  steered,  the 
Frost  algorithm  provides  excellent  control  over  the  response  in  the  chosen  look 
direction.  The  Frost-based  CBF  is  an  appropriate  choice  when  a  single  look 
direction  can  be  established  with  reasonable  accuracy. 

In  many  cases,  however,  it  is  not  possible  or  desirable  to  specify  a  single 
direction  of  arrival  for  desired  signals.  These  situations  arise,  for  example,  when 
communications  must  be  maintained  with  multiple  stations,  when  the  location  of 
a  station  is  imprecisely  known,  or  when  motion  of  the  receiving  platform  makes  it 
difficult  to  maintain  high  pointing  accuracies.  An  alternative  to  the  single,  hard- 
constrained  look  direction  is  needed  in  these  cases. 

There  are  two  requirements  that  must  be  met  in  order  to  realize  a  CBF  that 
provides  multiple  look  directions.  The  first  requirement  is  to  eliminate  signal 
energy  incident  from  the  various  look  directions.  It  is  difficult  and  expensive  to 
generate  wideband  nulls  in  multiple  directions.  On  the  other  hand,  it  is  relatively 


simple  to  create  multiple  narrowband  nulls  by  using  a  structure  such  as  a  nulling 
tree.  The  second  requirement  is  to  adaptively  eliminate  jammer  signals  while 
sustaining  useful  responses  in  the  look  directions.  Widrow’s  pilot-signal 
algorithm  provides  a  convenient  method  for  nulling  jammers  while  maintaining 
soft- constrained  responses  in  multiple  look  directions.  Taken  together,  the 
nulling-tree  preprocessor  and  the  pilot-signal  algorithm  form  the  building  blocks 
for  a  CBF  that  can  serve  multiple,  narrowband  signals. 

The  balance  of  this  appendix  describes  the  pilot-signal  CBF  and  draws  some 
comparisons  with  the  basic  pilot-signal  beamformer.  The  first  section  discusses 
the  nulling- tree  preprocessor  and  illustrates  the  use  of  the  nulling  tree  to  generate 
multiple  notches  and  broadened  notches.  The  next  section  then  explains  the 
CBF  structure,  which  essentially  consists  of  a  nulling-tree  preprocessor  plus  a 
pilot-signal  beamformer.  The  final  section  shows  that,  in  contrast  with  the  basic 
pilot-signal  ABF,  the  pilot-signal  CBF  is  capable  of  generating  an  unbiased 
solution  for  the  narrowband  case. 

A.  The  Nulling  Tree 

The  nulling  tree  was  described  by  Davies  [1]  and  is  also  termed  the  Davies 
cascade  or  Davies  null-steering  network.  As  has  been  noted  by  Steinberg  (2j,  the 
nulling  tree  does  not  necessarily  constitute  the  most  efficient  means  of  realizing 
multiple  nulls.  The  great  virtue  of  the  nulling  tree  lies  in  its  ability  to  provide 
independent  control  over  multiple  nulls.  In  that  sense  the  nulling  tree  is  well 
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matched  to  Widrow’s  pilot-signal  algorithm,  which  provides  independent  control 
through  the  separate  pilot  signal  established  for  each  look  direction. 

The  structure  of  the  nulling  tree  is  shown  in  Figure  D-l.  The  tree  consists 
of  a  number  of  stages  equal  to  the  number  of  nulls  to  be  generated.  Each  stage  is 
formed  from  a  set  of  identical  sections.  The  sections  operate  upon  two  outputs 
from  the  preceding  stage  and  perform  a  phase-shift- an d-add  process  that  nulls 
signals  arriving  from  a  selected  direction.  In  particular,  a  null  is  generated  in  the 
itk  stage  for  signals  arriving  from  an  angle  9  by  choosing  the  phase  shift  <j>{  given 
by 

<(>i  =  — D  -f  x  radians  ,  (D-l) 

where  d  is  the  interelement  distance  in  the  uniform  linear  array  and  X  is  the 
wavelength  at  the  frequency  for  which  nulling  is  desired. 

A  computer  program  was  written  to  design  nulling  trees  and  plot  the 
beampatterns  resulting  from  the  designs.  Various  test  cases  were  examined  to 
determine  the  characteristics  of  the  nulling-tree  preprocessor. 

Figure  D-2  shows  a  series  of  beampattern  pilots  that  result  from  adding 
successive  nulls  at  separate  arrival  angles.  The  initial  null  is  established  for 
signals  arriving  from  broadside;  the  second  stage  nulls  signals  from  +45*;  and 
the  third  stage  is  set  for  -30*  arrivals.  This  preprocessor  can  accommodate  an 
environment  in  which  there  are  multiple  signals  sources  at  well-defined  look 
directions. 


ling-tree  Beampattern  Plot  with  Nulls  at  -30  ,  0  ,  and  45 
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Another  type  of  requirement  is  considered  in  Figure  D-3.  Here  it  is  assumed 
that  only  one  emitter  is  of  interest,  but  the  look  direction  is  imprecisely  known. 
An  initial  null  is  established  in  the  estimated  look  direction,  and  successive  stages 
are  used  to  broaden  the  null  by  creating  additional  spatial  notches  on  each  side 
of  the  estimated  direction. 

B.  Structure  of  the  Pilot-Signal  CBF 

Figure  D-4  is  a  block  diagram  of  the  pilot-signal  CBF.  A  uniform  linear 
array  of  K  elements  provides  the  input  signals  to  the  CBF.  All  of  the  elements 
are  connected  to  the  nulling-tree  preprocessor,  which  consists  of  M  stages  that 
create  M  nulls  in  the  various  directions  from  which  desired  signals  may  arrive. 
The  K-M  outputs  from  the  nulling  tree  constitute  a  signal- free  environment  in 
which  adaptation  can  proceed  without  concern  about  cancellation  effects.  The 
slaved  beamformer  is  connected  directly  to  the  first  K-M  elements  of  the  array 
and  is  used  to  derive  a  useful  output  signal  for  the  receivers. 

The  pilot-signal  algorithm  is  used  to  provide  a  prespecified,  soft- constrained 
response  in  each  potential  look  direction.  A  separate  pilot  signal  is  generated  for 
each  look  direction  and,  after  appropriate  phasing,  is  added  to  the  ABF  inputs. 
The  pilots  are  also  filtered  and  summed  to  form  a  synthetic  “desired”  signal 
against  which  the  ABF  output  is  compared.  The  difference  between  the  summed 
pilots  and  the  ABF  output  is  the  error  signal  that  is  used  to  guide  adaptation. 


Figure  D-4.  Block  Diagram  of  the  Pilot-signal  CBF. 
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Weight  adaptation  can  be  accomplished  with  a  variety  of  different 
algorithms.  The  Widrow-Hoff  LMS  algorithm  [3]  was  initially  proposed  for  use  in 
a  digital  implementation  of  the  pilot-signal  ABF.  Applebaum  [4]  described  an 
analog  implementation  that  makes  use  of  correlation  loops. 

It  is  interesting  to  note  that  very  few  changes  are  needed  to  transform  the 
basic  pilot-signal  ABF  into  a  pilot-signal  CBF.  The  slaved  beamformer  is  already 
present  in  the  pilot-signal  scheme  and  therefore  need  not  be  introduced.  The 
only  new  components  are  the  null-forming  preprocessor  and  the  extra  antenna 
elements  entailed  by  the  nulling  process. 

C.  Bias  Removal  by  the  CBF 

Soon  after  the  development  of  the  pilot-signal  algorithm,  it  was  recognized 
that  the  solution  it  generates  in  the  presence  of  signal  is  biased  from  the  Wiener 
solution.  Griffiths  [5]  investigated  the  source  of  bias  in  some  detail  and  succeeded 
both  in  explaining  the  phenomenon  and  devising  an  alternative  beamformer  that 
avoided  the  problem.  Griffiths’  beamformer,  however,  is  subject  to  signal 
cancellation.  The  paragraphs  that  follow  will  demonstrate  that,  in  addition  to 
solving  the  signal-cancellation  problem,  the  CBF  solves  the  bias  problem,  thus 
eliminating  both  disadvantages  normally  encountered  with  pilot-signal 


beamformers. 
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Griffiths  investigated  the  bias  issue  in  the  context  of  a  pilot-signal  ABF  with 
a  single  look  direction  and,  consequently,  a  single  pilot  signal.  The  same 
approach  will  be  followed  here  because  the  use  of  multiple  pilot  signals 
complicates  matters  without  contributing  anything  of  substance. 

The  pilot  signal  in  Griffiths’  model  is  taken  to  be  some  multiple  0  of  the 


estimated  desired  signal  g(k),  i.e., 


rf(*)  =  0g(k)  . 


(D-2) 


After  phase  shifting  to  allow  for  propagation  effects,  the  pilot  is  added  to  the 
beamformer  inputs.  This  injected  signal  vector  is  identified  as  fiS(k),  where  S(k) 
is  the  estimated  desired-signal  vector  across  the  inputs. 

Weight  updating  is  accomplished  using  the  LMS  algorithm: 

IV(t+  1)  =  W(*)  +  2 nlPHk)  -  |Xr(t)  +  /Wr(*)|  W{k)}  {*(*)  +  /)5(t))  , 


(D-3) 


where 


X(k)  is  the  input  signal  vector, 

W(k)  is  the  beamformer  weight  vector,  and 
/i  is  the  adaptation  parameter. 

Griffiths  showed  that  the  limiting  value  of  the  expectation  of  the  weight 
vector  is  given  by 


where 


Rxx  4  E|X(Wr(*)l 
Rss  4  £|S(*)sr(*)l 
p,  4  £(*(*)$<*)]  ■ 

The  equation  for  differs  from  the  Wiener  solution 

W  =  RJ  P,  (D-5) 

even  if  the  pilot  happens  to  be  a  perfect  replica  of  the  desired  signal.  This  effect 
occurs  because  the  actual  incoming  signal  contributes  to  the  beamformer  inputs, 
and  the  combination  of  actual  signal  and  pilot  signal  enter  the  adaptation 
process.  Bias  in  the  pilot-signal  solution  is  undesirable  and  can  only  be 
eliminated  in  the  original  design  by  limiting  adaptation  to  periods  when  the 
desired  signal  is  known  to  be  absent.  This  is,  under  most  circumstances,  an 
unworkable  operational  constraint. 

In  the  pilot-signal  CBF,  adaptation  is  always  conducted  in  the  absence  of 
signal  due  to  the  action  of  the  preprocessor.  It  might  therefore  be  expected  that 
superior  performance  can  be  obtained.  This  is,  in  fact,  the  case. 

Some  further  examination  of  the  signal  environment  is  needed  to 
demonstrate  the  superiority  of  the  CBF  with  regard  to  bias.  The  signal  vector 
jf(k)  that  appears  at  the  beamformer  weights  is  taken  to  consist  of  two 
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components: 

*(*)  =  S(k)  +  N(k)  .  (D-6) 

The  vector  S(k)  Is  comprised  solely  of  desired-signal  components  arriving  via  the 
direct  path.  The  noise  vector  N(k)  is  the  sum  of  all  interfering  signals  (including 
desired-signal  multipath)  and  element  noise.  The  autocorrelation  matrix  for  X(fc) 
can  now  be  calculated: 

Rxx  4  £|X(*prT(*)l 

=  £{|S(*)  +  N(k)UST(k)  +  NT(k)}> 

=  £[S(t)Sr(t)|  +  £(S(t)Wr(t)l  +  £[iV(t)-Sr(i)]  +  £[JV(t)Wr(t) 1 
~  Rss  +  Rsn  +  vs  +  Rnn  P-7) 

Several  different  sets  of  assumptions  concerning  the  signal  environment  and 
the  CBF  realization  will  now  be  considered: 

Case  1:  Assume  that  signal  and  noise  processes  are  uncorrelated,  i.e.,  that 
multipath  and  repeater  jamming  are  not  present.  For  this  environment 

ft  XX  =  ftss  +  ftm  •  I0"®) 

Further  assume  perfect  implementation  of  the  CBF.  A  perfect  preprocessor 
excises  the  desired  signal  without  affecting  noise  components  and  generates  a 
modified  beamformer  input  given  by 


*'(*)  *  *(*) 


(D-9) 


The  autocorrelation  matrix  for  X'  (k)  is 


R^r  =  #, 


m 


(D-10) 


For  a  perfect  choice  of  pilot-signal  amplitude  and  statistics 


ft  =  1 
Rss  =  Rss 

p  —  p 
r»  ri 


(D-ll) 


The  limiting  value  KIW  for  Case  1  may  now  be  written  as 


nw  =  pEt+a. 


=  l#AW  +  #ssl-1  # r 


(D-12) 


This  last  equation  indicates  that,  unlike  perfect  realization  of  the  basic  pilot- 
signal  ABF,  perfect  realization  of  the  pilot-signal  CBF  yields  the  Wiener  solution. 


Case  2:  Repeat  the  assumptions  of  Case  1,  except  for  the  assumption  of  perfect 
preprocessing.  Instead  allow  perfect  removal  of  the  desired  signal,  but  assume 
incidental  modification  of  the  noise  field.  Under  these  conditions 

5=5  l-A  Rnn  +  #55 j-1  P f  ,  (D-13) 

where  the  A  matrix  describes  the  preprocessor  influence  on  the  noise  vector.  The 
right-hand  side  of  this  equation  is  not  necessarily  equal  to  the  Wiener  solution. 


Under  the  set  of  assumptions  discussed  in  Chapter  4  (i.e.,  narrowband 
interference,  sufficient  degrees  of  freedom,  high  jammer- to- noise  ratios  element 
spacing  less  than  half-wave),  the  right-hand  side  does  reduce  to  the  Wiener 
solution. 

Case  3:  Repeat  the  assumptions  of  Case  1,  except  for  the  assumption  that  the 
signal  and  noise  processes  are  uncorrelated.  When  correlation  is  allowed 

Rxx  —  Rss  +  Rsn  +  ^ns  +  Rm  •  (D-14) 

Preprocessing  now  produces  a  much  simpler  environment  within  the  adaptive 
beamformer: 

*'(*)  =  N(k) 

Rtx1  =  Rm  • 

The  limiting  value  of  the  weight  mean  is 

KJW  =  [ RpfN  +  Rss]-1  R 9  (D-16) 

In  other  words,  the  pilot-signal  CBF  treats  the  correlated  signal  as  though  it  were 
uncorrelated.  This  is  not  the  Wiener  solution,  but  it  is  the  desired  solution  in 
this  important  case. 
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APPENDIX  E 

SIMULATION  METHODS 

Computer  simulations  of  the  Frost  ABF  and  the  Frost-based  CBF  were  used 
to  illustrate  signal-cancellation  effects  and  to  demonstrate  the  performance 
contrasts  between  the  two  beamformers.  This  appendix  outlines  the  simulation 
methods  that  were  employed. 

A*  Simulation  Facility 

The  adaptive-antenna  simulations  for  this  research  were  performed  on  an  HP 
1000F  minicomputer  facility  at  Stanford  University.  The  simulation  software 
was,  for  the  most  part,  an  existing  Fortran  package  that  was  written  by  a  team 
of  students*  during  1979-1980  under  the  joint  sponsorship  of  the  U.  S.  Air 
Force’s  Rome  Air  Development  Center  and  the  Naval  Air  Development  Center. 
The  software  package  is  organized  as  a  set  of  programs  that  run  under  the 
control  of  a  driver  program  and  exchange  data  through  disk  files.  This 
arrangement  allows  multiple  simulation  functions  to  be  performed  in  parallel 
(albeit  on  a  time-shared  basis)  and  facilitates  research  by  allowing  time- 
consuming  adaptation  experiments  to  proceed  in  the  background  while 
experiment  preparation  and  data  analysis  are  being  performed  by  the  user.  The 
software  package  is  capable  of  simulating  several  different  beamformer  types  and 

*Priaeipal  eoatribiton  w*r»  R.  Gooch,  W.  Maaaiai,  F.  McCuthj,  W.  Ncwmoa,  aid  P.  Titckcaer. 
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is  sufficiently  flexible  to  support  most  of  the  basic  investigations  that  might  be 
required  for  a  specified  adaptive  antenna  and  signal/jammer  scenario. 

B.  Simulation  of  an  Adaptive  Antenna 

The  program  flow  for  adaptive  antenna  simulation  is  shown  in  Figure  E-l. 
The  initial  step  in  preparing  a  simulation  experiment  is  to  thoroughly  specify  the 
problem  that  is  to  be  considered.  The  driver  program  organizes  this  process  by 
prompting  for  the  inputs  that  are  needed  to  specify  the  following: 

A.  Antenna  configuration 

1.  Number  of  elements 

2.  Location  of  each  element 

B.  Beamformer  structure 

1.  Number  of  taps  per  beamformer  filter 

2.  Adaptive  algorithm 

3.  Constraints  (if  applicable) 

4.  Initial  weight  values 

C.  Signal/jammer  scenario 


1.  Number  of  signals /jammers 
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Figure  E-l .  Flowchart  for  Adaptive  Antenna  Simulation. 


.'A 


2.  Signal/jammer  arrival  angle 

3.  Signal/jammer  power 

4.  Signal/jammer  center  frequency  and  bandwidth 

D.  Adaptation  parameters 

1.  Number  of  adaptation  cycles 

2.  Adaptation  constant 

E.  Signal  (e.g.  input,  master  output,  slave  output)  to  be  stored  for  later 
analysis. 

All  of  these  data  are  organized  into  disk  files  and  thus  made  available  for  use  by 
various  programs  during  the  simulation  or  for  reuse  in  subsequent  runs. 

The  first  step  in  the  actual  simulation  is  generation  of  the  initial  samples  of 
the  input  signals.  This  can  be  accomplished  by  a  direct  computation  for 
sinusoidal  signals  since  amplitude  is  known  from  the  power  specification  and  the 
phase  at  any  tap  is  known  from  the  antenna/beamformer  specifications,  the 
signal  frequency,  and  the  initial  (or  prior-sample)  phase.  The  process  is  more 
difficult  for  wideband  signals.  These  signals  are  derived  by  first  generating  a 
white-noise  sequence  and  then  digitally  bandpass  filtering  to  obtain  the  desired 
center  frequency  and  bandwidth.  In  general,  the  sample  values  derived  in  this 
fashion  will  not  correspond  exactly  to  the  sample  values  needed  at  the  antenna 
elements  since  fractional-sample- interval  delays  are  often  encountered.  A  sine 


interpolation  is  used  to  correct  for  this  effect  and  provide  accurate  input  samples 
in  the  wideband  case.  The  aggregate  signal  sample  at  each  tap  is  just  the  sum  of 
the  various  signal  and  jammer  samples  derived  for  the  tap.  Separate  signal 
vectors  are  maintained  for  master  and  slave  beamformers. 

Once  the  signal  vectors  have  been  obtained  for  some  sample  instant,  it  is 
straightforward  to  calculate  the  output  signal  or  signals  by  forming  the  inner 
product  of  a  signal  vector  with  the  current  weight  vector.  In  the  case  of  the 
composite  beamformer,  both  the  master  beamformer  output  and  the  slaved 
beamformer  output  are  required. 

Weight  updating  is  accomplished  based  upon  the  current  signal  vector  and 
the  current  adaptive  array  output.  The  weight  recursion  for  the  Frost  ABF  was 
given  in  Eq.  (2-17). 

After  the  new  weight  vector  has  been  obtained,  the  next  adaptation  cycle 
can  begin.  A  check  is  made  at  this  point  to  determine  whether  the  prespecified 
number  of  cycles  has  been  completed.  If  so,  the  adaptation  routine  halts.  If  not, 
the  process  of  computing  input  samples,  outputs,  and  new  weight  values 
continues. 

As  the  simulation  progresses,  data  of  interest  are  saved  by  transcribing  to 
disk  files.  The  use  of  disk  files  allows  large  amounts  of  data  to  be  saved  and  thus 
supports  lengthy  adaptation  experiments  involving  tens  of  thousands  of 
adaptation  cycles.  Additionally,  the  disk  files  provide  convenient  interfaces 
between  the  simulation  program  and  various  data  analysis  programs. 
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C.  Plotting  of  Antenna  Responses 

Two  types  of  plots  are  typically  used  to  characterize  the  response  of  an 
adaptive  antenna  with  some  specified  set  of  beamformer  weights.  One  plot  shows 
the  response  to  a  unit-power  test  signal  of  specified  frequency  as  the  test-signal 
angle  of  arrival  is  varied  through,  say,  360*.  This  plot  is  usually  termed  a 
“beampattern”  and  is  useful  in  determining  the  relative  response  of  the  system  to 
jammers  and  desired  signals.  The  second  plot  shows  the  response  to  a  unit-power 
test  signal  from  a  specified  direction  as  the  test-signal  frequency  is  varied  over 
some  range,  generally  from  zero  to  half  the  sampling  frequency.  This  plot  is  a 
frequency-response  characterization  for  the  selected  direction  and  is  useful  in 
checking  the  effectiveness  of  jammer  nulling  and  confirming  the  maintenance  of 
look-direction  constraints.  Both  plot  types  are  used  most  frequently  to  study 
antenna  response  after  the  beamformer  weights  have  converged  (in  the  mean)  to 
the  optimal  solution  for  the  environment  under  consideration.  The  plots  may, 
however,  be  plotted  for  any  specified  set  of  weights.  A  series  of  plots  may  be 
generated  to  study,  for  example,  the  evolution  of  the  antenna  response  as 
adaptation  occurs  or  the  dynamics  of  the  antenna  response  after  convergence. 

A  flowchart  of  the  procedure  for  generating  and  plotting  antenna  responses 
is  given  in  Figure  E-2.  The  response  plots  are  specified  from  the  driver  program, 
which  then  starts  the  appropriate  program  for  either  a  beampattern  plot  or  a 
frequency-response  plot.  Either  program  must  first  read  the  driver-specified 
weight  file  and  plot-parameter  files  from  disk.  The  response  for  a  given  arrival 
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Figure  E-2.  Flowchart  for  Generation 

and  Plotting  of  Antenna  Responses. 


direction  and  frequency  is  then  found  by  forming  the  phasor  sum  of  the  weighted 
tap  voltages  and  computing  the  beamformer  output  power  from  the  amplitude  of 
the  sum.  Repeated  computations  of  this  type  give  the  response  at  enough 
directions  or  frequencies  to  accurately  characterize  the  aDtenna  response.  The 
response  is  then  plotted  for  study.  A  polar  plot  with  a  logarithmic  (dB)  radial 
scale  is  usually  preferred  for  beampattern  plots,  but  a  semilogarithmic  cartesian 
presentation  sometimes  serves  well.  A  cartesian  plot  is  used  for  the  frequency- 
response  data,  with  frequency  as  the  abscissa  and  power  (in  dB)  as  the  ordinate. 

D.  Plotting  of  Input/Oatput  Signals 

Plots  of  input  and  output  signals  and  their  spectra  are  important  in  the 
study  of  adaptive  antenna  performance.  These  plots  are  especially  valuable  in 
investigations  of  signal  cancellation.  Cancellation  effects  that  are  difficult  to 
recognize  in  beampattern  or  frequency-response  plots  can  be  easily  detected  by 
comparing  input  and  output  waveforms  or  spectra. 

Input  and  output  data  are  gathered  during  the  course  of  simulation 
experiments  and  written  to  disk  files  for  later  analysis  and  display.  A  number  of 
different  signals  (e.g.  the  desired  signal,  jammer  signals,  master  beamformer 
output,  slaved  beamformer  output)  may  be  collected  to  thoroughly  characterize 
beamformer  behavior  during  various  simulations.  Plots  of  these  data  allow 
comparisons  to  be  drawn  between  different  beamformers,  different  operating 
conditions,  etc. 
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Figure  E-3  diagrams  the  event  Sow  for  production  of  input/output  plots.  If 
spectral  plots  are  elected,  a  DFT  routine  is  entered  to  generate  power  spectra  for 
the  selected  signals.  The  DFT  routine  allows  selection  of  transform  size  and 
supports  the  computation  of  an  average  power  spectrum  over  some  specified 
range  of  data  samples.  The  power  spectra  are  written  to  disk  upon  completion  of 
this  computation.  The  plot  program  allows  several  signals  or  spectra  to  be 
presented  simultaneously  for  comparison.  Waveform  or  spectral  data  are  read 
from  specified  disk  files  and  plotted  on  separate  axes  with  a  common  scale  for 
easy  comparison. 
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Figure  E-3.  Flowchart  for  Generating  and 
Plotting  Input/Output  Signals 
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